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GALACTIC AND EXTRAGALACTIC STUDIES, XX. EXTENT 
OF THK GALAXY IN THE ANTICENTER OCTANT 


By HARLOW SHAPLEY AND ANN B. HEARN 
HARVARD COLLEGE OBSERVATORY, CAMBRIDGE, MASS 
Communicated July 31, 1952 


Further application of the galaxy-cepheid method of estimating sidereal 
distances in and near the Milky Way has made possible this report on the 
extent of the galactic system in the anticenter region. The discovery and 
study of the variables and galaxies in the selected area, which encom 
passes five thousand square degrees, bounded by 


8 (galactic latitude) +40° to —40 
\ (galactic longitude) 110° to ISO 


is far from complete, but in view of the increasing interest in the total di 
mensions of galaxies the results now on hand are summarized below. 

1. Eighty-four cluster-type cepheids have now been found in the area 
and when possible their distances in kiloparsees have been computed from 
the usual relation 


d - 10° 2im Am VW 2 


where my) is the observed median magnitude derived from the light curve 
and Am is the absorption in magnitudes between observer and star. By 
adopting the usual value of the median absolute photographic magnitude, 
M = 0, and using the corrected median apparent magnitude, rm, = rte 

Am, we have log d = 0.2m, — 2. The distance of a variable star from the 


galactic plane, its distance from the observer projected on the center-anti 


center line in the galactic plane, and the distance from the center of the 
galaxy to the star are given, respectively, by 
z=dsn£B 


w = d cos B cos (A 145 
d, = [22 + s? + (d cos B)? + 2sw 


where 8 and X are the galactic coordinates, the longitude of the galactic 
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center 1s taken as 325”, and s is the distance from the observer to the galac- 
tic center. For the following computations we shall use both s = 8.3 
kpe. (d,(1) in table 1) and s = 9.0 kpe. (d,(2) in table 1), the latter value 
based on an unpublished study by Mr. Thomas A. Matthews of the Har- 
vard Observatory. ! 

2. Toestimate the absorption, 4m, we have assumed that the true popu 
lation of galaxies in the anticenter octant, when freed of interstellar ab- 
sorption, is comparable with the average population for the high latitude 
sky. To the extent that this assumption fails, through significant irregu- 
larities in the distribution of faint galaxies, to that extent also are the in- 
dividual distances uncertain. The general order of the distances, we be- 
lieve, is dependably estimated and therefore a true indication of the galactic 
extent can be derived from this material. We have, of course, recognized 
the probability that our standards of average galaxy population, which are 
in the galactic polar zones, are also not free of absorption and therefore for 
all of the anticenter octant a correction of 0.25 magnitude is added to what 
ever absorption is shown by the relative deficiency in the population of 
faint galaxies along the Milky Way. This correction has decreased all 
computed distances by eleven per cent. 

4. During the past two years the coverage of the anticenter region with 
long exposure plates for the detection of faint galaxies has been brought to 


seventy-five per cent completion. Two Agassiz station telescopes, the 16 


inch Metealf and the 24-33-inch Jewett-Schmidt, have contributed 
most of the new long exposure photographs. The shorter exposure plates 
for the discovery and study of variable stars are mostly made with the Met- 
calf instrument. 

From the examination of the long exposure plates we derive the absorp 
tion through the relation connecting galaxy numbers, V, with apparent 
magnitudes 


log N = b(m) — Am — m) 


With b = 0.6, m, = 105.2, as the paramciers for uniformly dense population 
of galaxies throughout the regions here explored, 


Am = my ~— 15.2 — 1.67 log N 


where for each plate my is the observed faintest magnitude for galaxy im 
ages. From this relation we have computed the absorption for the deter- 
mination of m, and the distances for the 52 stars in table 1. (Four stars of 
VSF 211 and 309 are slightly outside the bounds of the octant.) 

!. For $2 additional cluster variable stars the galaxies in the surround- 
ing fields are invisible, or so few, faint, and uncertain that no attempt has 
been made to evaluate the absorption and distances. For 11 stars in 
table 1, where Am is one magnitude or greater, the distances are neces- 
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TABLE 1 


Firry-Two VARIABLES IN THE ANTICENTER Oct: 


NAME -E RIOD 


RX Eri 5872 
X Ari 6511 
BC Eri 2639 
RR Gem 3973 
TT Cne 634 
SZ Gem 5011 
AH Tau 1663 
$8. 1933" 
SS Tau 3700 
TZ Aur 3917 
SS Cne 3673 
RV Ari O852 
RW Ari 2614 
rY Cam 6701 
44.1929" 
HV 7668" 
RZ Cam 
XY Eri 
SX Cne 
HV 10389 
HV 7649" 
GP Aur 
HV 7661" 
YZ Tau 4115 
CR Aur 5914 
BE Eri 5795 
AL Eri 6569 
CV Tau 
AC Eri 1821 
BN Eri 
BG Eri 6599 
BF tri 
AI Tau 5686 
HV 10416' 
WY Eri 5069 
HV 7665 
AG Eri 5125 
XX Eri 7102 
CK ‘Tau 
AE Eri O514 
AB Eri 5740 41 
GO Ori 21 
HV 6284 
Al Eri 6195 30 
AP Eri 4211 37 
AD Eri 0 4699 10 
HV 10402 30 
HV 7655" 16.5 
BO Eri 24 
HV 10396 23 
VZ Eri 35 43 
FV Ori 52: 24 


Vsti wewenwean cos 


ao 


H 
>.2 
1.4 
5.8 


Nae DeS DN 


* Not certainly classified as a cluster-type cepheid 
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sarily uncertain. For the brightest cluster variables, and this applies also 
to the 65 classical cepheids which we do not use because of their small dis 
tances or low latitudes, we cannot be sure, in the absence of measures of 
color, whether the variables are actually beyond the absorbing material 
or init, or in front. We may have overcorrected. But for all the regions 
with 6B > +10°, where faint galaxies show, we can assume safely that the 
variables lie far beyond the absorption and the full value of Am should be 
applied. 

The median magnitudes of the variable stars and the plate limits for the 
galaxy survey are based on various magnitude standards, including selected 
area sequences and, for the higher latitude fields, sequences set up through 
star counts. They are of course subject to the usual uncertainties of 
photographic photometry, but the errors probably average out; in any, 
case, because of the inescapable uncertainties in the correction for space 
absorption, the errors may be neglected for this qualitative investigation. 

5. The variables are listed in table | in order of increasing distance 
of the star from the center of the galaxy. It is at once evident that these 
stars are members of our galactic system's stellar corona and not members 
of the rotating galactic discoid. 

The low latitude obscuring clouds have so far prevented us from measur- 
ing the total extent of the discoid in the anticenter region. We have as yet 
little evidence that the discoid can be traced more than a few kiloparsecs 
in the Taurus-Auriga direction, although some of the fainter classical ceph- 
eids in Monoceros as well as in Cygnus suggest that such discoid-inhabiting 


objects are seen at great distances from us and from the galactic center. 


6. Although we have not as yet measured any discoid members at large 
distances in the anticenter region, we have in this report shown the great 
extent of the stellar corona in the anticenter octant. It extends consider- 
ably beyond the radius of 15 kpe that we have in recent years commonly ac 
cepted for the galactic system's radius. For 12 stars the distance from the 
galactic center (d,(1)) is greater than 50,000 light years; or for 15 stars 
(207), if the larger value of s is adopted. 

If the cluster variables provide a fair indication of population density we 
find from the evidence at hand a very low stellar density at the farthest 
distances here explored. The contribution to the total galactic mass from 
the outer corona must be quite negligible. A similar conspicuous decrease 
of population density with distance from the galactic plane, particularly in 
the direction of the south galactic pole, was reported in an earlier com- 
munication.2. We are now making a comparable study of the cluster 
variables in the north galactic polar zone and have also initiated a program 
of exploring the corona at intermediate latitudes. 

The Office of Naval Research is currently assisting in our study of the 
anticenter section of the Milky Way. 
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Summary. —The extent of the corona of stars that surrounds the galactic 
system has been examined on the basis of current studies of galaxies 
and variable stars in the anticenter region of the Milky Way. The vari 
ables of the cepheid family found in galactic latitudes of less than ten degrees 
can rarely be used to measure the extent of the Milky Way because the 
interstellar absorption is heavy; but on the borders of the Milky Way, 
where the galaxy counts indicate the degree of transparency, a considerable 
number of cluster type variables have shown that the coronal population 
extends to anticenter distances equal to the distance separating us trom the 
center of the Milky Way. The distances from the center of the galaxy to 
a fourth of these outlying variable stars exceed fifteen kiloparsees; similar 
distances were found previously for a few variables in the galactic polar 
caps. The greatest distance exceeds twenty kiloparsecs. The computed 
distances above and below the galactic plane (z in table |) show that few if 
any of these cluster variables are members of the galactic discoid, 

His mean value of s is 9.0 + 0.6 kpc 
* These PROCEEDINGS, 25, 423-428 (1939); Harvard Reprint 173 


NATURAL AND ACQUIRED RESISTANCE OF ESCHERICHIA 
COLI TO NEOMYCIN* 


By WARREN P. IVERSON AND SELMAN A. WAKSMAN 


DEPARTMENT OF MICROBIOLOGY, NEW JERSEY AGRICULTURAL EXPERIMENT STATION, 
RUTGERS UNIVERSITY 


Communicated August 11, 1952 


Introductory... The problem of development of resistance of bacteria and 
other organisms to the action of bacteriostatic and bactericidal agents in 
general and to chemotherapeutic substances in particular has become in 
creasingly serious with the extensive use of antibiotics in the treatment of 
infectious diseases of men and ammals. Different compounds exhibit 
marked differences in the nature and speed of development of resistance, 


thus pointing to differences in their mode of action upon bacteria and other 
microorganisms. Some antibiotics, like streptomycin and more especially 


grisein, favor rapid development of resistance among sensitive organisms. 
Others, like penicillin, favor development of resistance among various bac 
teria only upon continued use. Moreover, penicillin leads to the develop 
ment of resistance only among certain bacteria and not at all or to a lesser 
degree among other organisms; the rate with which such resistance de 
velops and the degree of resistance attained are not so great as for strep 
tomycin. Aureomycin and terramycin do not allow the rapid develop- 
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ment of resistance among sensitive bacteria, although these antibiotics 
may be reduced in potency in certain media. Still other antibiotics, like 
aspergillic acid and polymyxin, are said not to allow any development of 
resistance. 

Among the antibiotics which have received particular consideration both 
because of their importance as chemotherapeutic agents and because of 
the type of resistance developed by them, streptomycin occupies a promi 
nent place. This antibiotic favors not only rapid development by various 
bacteria of a certain limited degree of resistance, from less than | unit or 
microgram per milliliter to 20 wg./ml., but also rapid attainment of abso 
lute resistance, to that of more than 5000 wg./ml., thus frequently render 
ing this antibiotic almost valueless upon prolonged administration. 

Since the first announcement, in 1944, of the isolation of streptomycin, 
numerous organisins, mostly actinomycetes, have been examined for their 
antibiotic-producing capacity in an effort to find agents that would possess 
antimicrobial activities similar to those of streptomycin, especially against 
gram-negative and tuberculosis bacteria, but without favoring the rapid 
development of resistance as a whole and of absolute resistance in particu 
lar. Neomyein appeared to be one such agent. In the first report of the 
isolation of this antibiotic,’ 1t was emphasized that neomycin is effective 
against both streptomycin-sensitive and streptomycin-resistant bacteria. 
Out of a population of 246 million cells, S units of neomycin brought about 
complete inhibition of growth, 6 units permitted the development of only 2 
colonies, and 4 units allowed the growth of LOS colonies, after 5 days’ incu 
bation; corresponding concentrations of streptomycin yielded 70, 130 and 
640 colonies, respectively 

In further studies® of the comparative bacteriostatic and bacteriolytic 
properties of streptomycin and neomycin, the latter was found to possess a 
greater bactericidal effect upon /¢scherichia coli than the former. Com 
binations of the two antibiotics showed an additive and possibly even a 
synergistic activity. Streptomycin-resistant strains of /. coli showed no 
increase in resistance to neomycin, but neomycin-resistant strains showed 
a slight increase of resistance to streptomycin. Similar observations were 
made by Demeree and Demerec? for /:. coli and by Gezon and Fasan‘ for 7 
strains of group AB hemolytic streptococci. They found that neomyein 
shows a resistance pattern similar to that of penicillin. Yegian and Van 
derlinde® reported, however, that although neomycin is just as effective 
against streptomycin-sensitive tubercle bacilli, and vice versa, the incidence 
of resistance is similar for both antibiotics. 

In view of the growing importance of neomycin as a chemotherapeutic 


agent in the treatment of skin infections and by oral administration, a 


iurther study of the problem of resistance to this antibiotic was felt justi 


fied. 





Vor. 38, 1952 BACTERIOLOGY: IVERSON AND WAKSMAN 845 


Experimental. — Sensitivity of Different Strains of k. coli to Neomycin: 
Nutrient agar plates containing different concentrations of neomycin were 
streaked with aqueous suspensions of several strains of /:. coli grown for 
20 hrs., at 28°C., on nutrient agar slants. After incubation for 2 days, the 
plates were examined for bacterial growth. The results presented in table | 
show that, with the exception of one strain, all the cultures of /:. coli were 
completely inhibited by a concentration of 5 u/ml. of neomycin. The 
strain (/:. colt 57) that was found to be resistant to neomycin proved to be 


TABLE 1} 


SENSITIVITY OF VaRIOUS FE. coli STRAINS TO NEOMYCIN 
NEOMYCIN IN PLATES, U) ML. OF AGAR*® 
STRAIN O95 l 


colt 56 y | 


coli 5s 


1 

cola 57 } 
: | 

| 

| 


E 

I 

F 

E. coli 54 

E c olt 55 

'E. colt 52 

FE. coli 52SR ; 2 
FE. coli 52SD 0 


excellent growth; 38 = good growth; 2 = fair growth; | = poor growth 


> 


7 = trace of growth in the form of isolated colonies. SR streptomycin-resistant 
SD = streptomycin-dependent 


* ” 


lo convert units of neomycin to micrograms, multiply by 3.8 


rABLE 2 
COMPARATIVE INHIBITORY ACTION OF NEOMYCIN AND STREPTOMYCIN ON EF. colt 
Survivors per ml. cell suspension 


ANTIBIOTIC 
| MI AGAR STREPTOMYCIN NEOMYCIN 


7.400 000 1700 
H40 108 
130 2 
70 <1 
1 ml. suspension contained 246 & 10° cells, 24 hrs. old Figures are the average 


colony counts of 4-6 plates 


resistant also to streptomycin. The streptomycin-resistant strain (/. 
coli 52SR) was only slightly more resistant to neomycin than was the 
streptomycin-sensitive strain. The streptomycin-dependent strain made 
no growth in the presence of neomycin.® 

Presence of Neomycin-Resistant Cells in a Normal Culture of E. colt: 
A suspension of washed cells from a 24-hr.-old nutrient agar slant of /. coli 
gave 246 X 10° viable cells per milliliter. The ceil suspension was diluted 
with sterile water and plated out in nutrient agar containing either strep 
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tomycin or neomycin in concentrations of 2, 4, 6 and 8 units per ml. The 
plates were incubated for 5 days at 28°C., and the colonies developing on 
the plates counted (table 2), Streptomycin allowed the development of a 


great taany more colonies in a given suspension of F. coli cells than did 
neomycin. For example, 2 u/ml. of neomycin brought about the inhibi- 
tion of about 4300 times as many cells as did 2 u/ml. of streptomycin. 
With higher concentrations of the antibiotics this ratio was reduced, but 
the difference was still significant. 


TABLE 4 


NEOMYCIN RESISTANCE OF 5 CULTURES OF F. coli ISOLATED FROM COLONIES GROWING 
IN THE PRESENCE OF 6 U/ML. OF NEOMYCIN 


NEOMYCIN 
! Mt 4 - 


COLONY NO 
) 3 


Growth 


rABLE 4 
DEVELOPMENT OF NEOMYCIN RESISTANCE IN A CULTURE OF EF. coli GROWING IN NUTRI- 
ENT Brotu with 0.4 U/ML OF NREOMYCIN 


POTAL NUMBER OF CELLS RESISTANT TO 


INCUBATION 
VIABLE CELLS/ MI 5 U/ML. OF NKOMYCIN 


HRS 
0 2,280,000 <1 

j $21 3 
s 91 8 
12 247 28 
2A 19,500 1,640 
18 63,000 , 000 420 , O00 
72 210,000,000 1,000 ,O000 
06 241,000 ,O000 20 ,000 ,000 
144 293 ,000 ,000 760 ,000 
219 220 000 ,000 210,000 


In another experiment, a suspension of -/:. coli containing 22 billion cells 
per milliliter was plated out in an agar medium containing 10 u/ml. of 
neomycin and a suspension of 44 billion cells in a medium containing 25 
u/ml. of neomycin. After 9 days’ incubation, no colonies of neomycin- 
resistant or -dependent organisms appeared on any of the plates. 

Development of Neomycin-Resistant Cells During the Growth of E. coli: 
An effort was made to isolate colonies more resistant to neomycin than was 
the bacterial population as a whole. An £. co/i suspension containing 246 
million viable cells per milliliter was plated out on nutrient agar having a 
final concentration of 6 u/ml. of neomycin. <A total of eight colonies ap- 
peared on the plates. When transferred to nutrient agar, only five of the 
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colonies gave growth. The other colonies failed to produce any growth on 
transfer to fresh agar medium with or without neomycin. The cultures 
obtained from the five colonies that produced growth were tested for their 
sensitivity to neomycin, by means of the agar streak method (table 3). 
These cultures were found to be only slightly more resistant to neomycin 
than was the original culture. Microscopic examination of the cultures 
growing on neomycin-containing agar showed the cells to be elongated, 
being two to three times the length of the cells on neomycin-free agar. 
The growth of the cultures in normal broth was similar to that of the parent 
strain, both in size and in general appearance. 

Next, a study was made of the development of resistance to neomycin 
in a culture of /. coli growing in nutrient broth containing this antibiotic. 
Three hundred ml. nutrient broth with 0.4 u/ ml. of neomycin, an amount 
insufficient for complete growth inhibition, was inoculated with a fresh 
culture of /. coli. The inoculated broth culture was incubated at 28°C., 


rABLE 5 
DEVELOPMENT OF STREPTOMYCIN RESISTANCE IN A CULTURE OF FE. colt GROWING IN 
BROTH WITH 5 U/ML. OF STREPTOMYCIN 


INCUBATION TOTAL NUMBER OF CELLS RESISTANT TO STREPTOMYCIN, wG MI AGAR 
HRS VIABLE CELLS/ML 4 8 12 


0 1,980,000 2%; , 
: 53 108 4 
S 900 347 , 
19 $000 1.357 149 
18 46,000,000 2,310,000 180.000 180 
144 239 ,000 ,000 239 ,000 ,000 239 000.000 239 000,000 
189 606 , O00 O00 634,000 000 623 ,Q00 000 569,000,000 


and aliquots were taken at frequent intervals and plated out on ordinary 
nutrient agar as well as on agar containing 5 u/ml. of neomycin (table 4). 
A similar study was made of the development of cells resistant to strepto 
mycin in the presence of 5 u/ml. of this antibiotic in the broth (table 5). 
The addition of 0.4 u/ml. of neomyein was about as effective as 5 u/ml. 
of streptomycin in reducing the number of viable cells during the first 4 
hrs. of growth of a fresh culture of /. coli. After 24 hrs.’ incubation, the 
number of viable cells greatly increased. The increase in resistance to 


streptomycin was far more rapid, however, than that to neomycin. Al 
though only | to 10 per cent of the cells became resistant to neomyein (5 u 
inl.), nearly all became resistant to streptomycin (12 u/ml.). 

Reversion to Sensitivity of Neomycin-Resistant Strains: One of the neo- 
mycin-resistant strains (No. 5) isolated in one of the previous experiments 
was selected for further study. When grown on nutrient agar, this strain 
produced two types of colonies that were almost indistinguishable except 


for their size. Type A colony was similar to that of the parent strain, the 
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individual cells being normal in appearance. Type B colony was smaller 
than that of type A, the cells being irregular in length, some forming long 
filaments with fusiform bodies. Type B appeared to be unstable when 
grown on neomycin-free agar, giving rise to a mixture of A and B types. 
The parent strain of /:. coli, the colony isolate, and the two substrains were 
tested for their sensitivity to neomycin, as shown in table 6. 

An aqueous suspension of the colony isolate No. 5 (644 & 10° viable 
cells/ml.) was plated out on nutrient agar containing 15, 20 and 25 u/ml. 
neomycin. After incubation for S days, a total of 59 colomes appeared 

FABLE 6 
NEOMYCIN SENSITIVITY OF TWO SUBSTRAINS DERIVED FROM A COLONY ISOLATE OF &, 
coli ORIGINALLY RESISTANT TO 6 U/ML. OF NEOMYCIN 
CULTURE 


SUBCUL SUBCUL 
TURE A TURE B 


NKOMYCIN 
usMi NO.5 


Growth 


rABLE 7 


NEOMYCIN RESISTANCE OF S SUBSTRAINS ISOLATED FROM COLONIES GROWING IN 
PRESENCE OF 15 U/ML. OF NEOMYCIN 


CULTURE NO 


NKOMYCIN 
5d Se 


u/MI Vi 5b 
Growth 

) 
15 


20 


4 
l 
l 
l 


25 


30 


—-w Ww Ww + 


0 


Many large isolated colonies 


on the plate containing 15 u/ml. neomycin, but none on the plates contain 
ing higher concentrations of this antibiotic. Fourteen of the colonies were 
isolated on slents. Only eight of these grew. They were tested for their 
sensitivity to neomycin (table 7)... All the substrains grew on the neomy- 
cin concentration from which they were originally isolated; all except two 
were completely inhibited by 30 u/ml. of neomycin. Of the two more re- 
sistant substrains that have developed, No. 5d was selected for further 
study. 

One-ml. portions of a suspension of No. 5d (563 X% 10° viable cells/ml. ) 
were plated out on nutrient agar containing 30, 35, 40, 45 and 50 units of 
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neomyein per milliliter. After 14 days’ incubation, no colomes appeared 
on any of the plates. When retested, the strain was found to have de 
creased in resistance, so that 20 u/ml. of neomycin was sufficient for its 
complete inhibition. 

Summary.—-Eight cultures of /. coli, including two streptomycin-re 
sistant strains, were tested by the agar-streak method for their sensitivity 
to neomycin. All cultures, except one, were found to be partly inhibited by 


0.5 u/ml. of neomycin and completely inhibited by 5 u-ml. One culture 
was naturally resistant to at least 10 u/ml. of neomycin. Large inocula 
required larger concentrations of neomycin for growth inhibition than did 


smaller inocula. In a natural population of /:. col: cells found in a given 
culture, a small number of cells were somewhat more resistant to neomycin 
than was the population as a whole. A much smaller number of cells were 
resistant to neomycin than to streptomycin. In the process of growth of 
FE. coli, the number of cells less sensitive to neomycin increased during the 
logarithmic phase and tended to disappear ducing the stationary phase of 
growth. No highly neomycin-resistant or neomycin-dependent cells were 
found in natural populations of /. coli. Cultures of /. coli that had de 
veloped resistance to neomycin showed a tendency toward reversion to 
sensitivity. Growth of slightly neomycin-resistant strains of /. colt in 
the presence of non-inhibitory concentrations of neomycin was very poor 
and was accompanied by marked changes in the morphology of the cells. 
* Paper of the Journal Series. New Jersey Agricultural Experiment Station, Rutgers 
‘niversity ~The State University of New Jersey, Department of Microbiology 
1 Waksman, S. A., and Lechevalier, H. A., Science, 109, 305-307 (1949) 
? Weiss, D., and Waksman, S. A., Proc. Nati. Acap. Scr., 36, 293-300 (1950) 
§ Demerec, M., and Demerec, R., Proc. SOth Ann. Meet. Soc. Amer. Bact., (1950), p 
101 
‘ Gezon, H. M., and Fasan, D. M., Science, 114, 422 423 (1951) 
> Yegian, D., and Vanderlinde, R. J., Am. Rev. Tuberc., 61, 488-507 (1950) 
® To make the units of measurement comparable for both neomycin and streptomycin, 
the unit is used as a standard. It is to be remembered that a unit of streptomycin is 


equal to | wg. and a unit of neomycin to 3.3 yg 





BIOCHEMISTRY: A. J. HODGE Proc. N. A. S 


A NEW TYPE OF PERIODIC STRUCTURE OBTAINED BY RE 
CONSTITUTION OF PARAMYOSIN FROM ACID SOLUTIONS* 


By A. J. Hopcet 


DEPARTMENT OF BroLoGy, MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


Communicated by F. O. Schmitt, August 1, 1952 


The type of periodic distribution of density characteristic of naturally 
occurring paramyosin fibrils (native paramyosin) obtained from molluscan 
adductor muscles has been described by Hall, Jakus and Schmitt.! Small 
angle x-ray diffraction work,’ using adductor muscles of the clam, Venus 
mercenaria, showed that the native paramyosin fibril has small angle order 
in only one plane. The periodicity in the dry state is described by a non 
primitive cell with elements do = 261 A., bo) = 720 A. (fiber axis), included 


angle 90°35’ and with five equivalent positions per cell. In the wet 


state, the structure is similar, but with cell elements a9 = 310 A., bp = 725A. 
and included angle 94°. The type of lattice deduced from the x-ray dif 
fraction data can be directly visualized in electron micrographs of fibrils 
stained with phosphotungstic acid (PTA), the structure consisting of fine 
dense bands spaced 145 A. apart, together with a network of dense spots 
superimposed on the bands in such a manner that the true repeat period is 
5 ™& 145 = 725 A. Thus, the native paramyosin fibril can be regarded as 
a well-characterized entity, easily recognizable in the electron microscope. 
Extensive electron microscopic investigation of preparations from Venus 
adductor muscles has failed to disclose the presence of any type of periodi 
cally striated fibril other than that described above.’ 

It will be shown in the present paper that it is possible under certain 
conditions to cause at least one of the constituents of native paramyosin 
(which itself may be a complex) to undergo a reconstitution process from 
solution, resulting in the formation of protein fibrils having a type of re 
peating structure distinetly different from that of the native paramyosin. 

When adductor muscles from Venus (particularly the white portions 
thereof) are macerated by means of a Waring Blendor in dilute salt solution, 
(0.2.0.3 47 KCI, pH 5-6), a suspension of long needle-like fibrils is obtained 
which, however, contains considerable amounts of dissolved proteins and 
cellular debris. Such suspensions are easily purified by differential cen 
trifugation, and after repeated centrifugation, it is possible to obtain sus 
pensions containing very little other than typical native paramyosin fi 
brils. These fibrils are insoluble at low tonic strength in the pH range 5-7. 
On the addition of dilute acid (usually acetic acid, because of its buffering 
capacity), there is an immediate clearing of the opaque suspensions as the 
protein goes into solution, Such solutions have ultra-violet absorption 
spectra characteristic of proteins, no unusual bands being observable. By 
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high-speed centrifugation or filtration through sintered-glass filters, it is 
possible to obtain clear solutions apparently free of fibrillar material within 
the limits of resolution of the electron microscope. A translucent gel of 
very fine filaments is obtained on prolonged dialysis of such solutions 
against large volumes of distilled water. These filaments are too thin to 
determine whether or not they possess periodic structure. However, if 
the pH is maintained in the region of 3.5 and the ionic strength of the 
dialysis fluid is increased to about 0.35, a fibrous precipitate of hitherto un- 
observed type is obtained. The periodic structure characteristic of this 
new type of fibril is illustrated in figures | and 2. The structure is sym 
metrical with respect to the disposition of intraperiod banding, is devoid of 


FIGURE 1 
PLS fibrils obtained by dialysis of a typical paramyosin solution against 0.4 7 KCI 
at pH about 3.5. Stained with 1% pbosphotungstic acid (PTA). Magnification 
71 000 X. 


densely staining spots such as are seen in the native fibril, and the fiber 
axis repeat period is in the vicinity of 1400 A. In conformity with the 


nomenclature adopted for an analogous situation in which a new long-spac 
ing structure was obtained from collagen preparations,' the new type of 
periodic structure will be designated as paramyosin long-spacing (PLS). 
The distribution of density within the repeat period of PLS bears a striking 
resemblance to that present in the sarcomere of vertebrate striated muscle, 
the PLS analogs of the Z, /, A, J and // bands all being readily discernible 
in electron micrographs (Fig. 1). Fibrils of PLS formed at pH ca. 3.5 
tend to be characteristically tactoidal and very thin. From shadow lengths 
the thicknesses were estimated to lie in the range 40-100 A. The periodic 
distribution of thickness observed in shadowed specimens coincides with 
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the density distribution in specimens stained with PTA (compare Figs. 1 
and 2). 

PLS fibrils can be obtained more consistently and in greater abundance 
by a citrate extraction procedure, which makes it unnecessary to isolate 
the native paramyosin fibrils in order to obtain suitable solutions. Whole 
minced adductor muscle (preferably white) is extracted overnight in the 
cold with M /6 potassium citrate buffer at pH 3.4. The residue is recovered 
by filtration or centrifugation, washed with the buffer and finally extracted 


FIGURE 2 
PLS fibrils from a preparation similar to that of figure 1. Stained with 1°, PTA 
and shadowed with chromium. Magnification 33,000 x. 


with dilute acetic acid (usually 0.4%, by volume) for several hours. The 
resulting viscous solution is cleared of fibrillar material by centrifugation or 
filtration as before. The PLS precipitates obtained from such solutions 
are generally more gel-like and less flocculent than those obtained from 
solutions of isolated native paramyosin fibrils and are thus more suitable 
for the preparation of specimens for examination in the electron micro- 
scope. 


The physical state of the precipitate depends markedly on the pH and 


ionic strength of the dialysis medium as well as on the method of prepara- 
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tion of the paramyosin solutions. The ionic strength required to produce a 
precipitate decreases with increasing pH, the precipitate being more gel- 
like and the constituent fibrils thinner. This tendency is sufficiently pro- 
nounced at pH above about 4 that it is difficult to observe the periodic 
structure characteristic of PLS. It has not yet been established whether 
this apparent lack of periodic structure is due to (a) the obvious difficulty of 
observing periodic structure in very fine fibrils or filaments, (b) the necessity 
for the formation of periodic structure of maintaining the dialysis condi- 
tions within certain limits or (c) the operation of both the above factors. 
As has been pointed out already, a gel of fine apparently structureless fila- 
ments is obtained on dialyzing the paramyosin solutions against distilled 
water. However, all attempts thus far to reconstitute fibrils having the 
native paramyosin type of periodic structure have failed. Nevertheless 
the possibility exists that such fibrils may be obtained from the solutions 
under certain conditions not yet attained. Should this not be possible, the 
evidence would suggest that the protein undergoes some physical change 
such that, having once been exposed to acid conditions, it is incapable of 
forming the native type of structure. 

In order to elucidate the properties of the paramyosin system, particu- 
larly with respect to its ability to give rise to a modified periodic structure, 
a physical chemical investigation has been undertaken. Preliminary data 
have been obtained for ultracentrifugal, diffusion, electrophoretic, viscosity 
and light scattering analysis, and work is continuing along these lines. 
Solutions capable of yielding PLS, whether derived from the citrate method 
or from isolated native paramyosin fibrils, have been investigated mainly 
with respect to variation in ionic strength. 

All methods of analysis so far utilized have indicated the presence, in 
the solutions, of highly asymmetric particles of relatively high particle 
weight. The solutions have high values for the viscosity increment, the 
experimental results indicating asymmetry ratios of the order of 50:1 in 
low ionic strength media. That the observed high viscosities under these 
conditions are not due to the electroviscous effect is suggested by the fact 
that the viscosity of the solutions at very low ionic strength is essentially 


independent of temperature between 0 and 40°C. Hence the observed 


high viscosity increment under these conditions is probably not due to 
particle interaction. In media of ionic strength 0.1 or higher, the solutions 
are at least SO™ homogeneous in the ultracentrifuge, the pattern consisting 
of a large well-defined peak and a much smaller second peak sedimenting 
more rapidly. The latter component does not seem to be necessary for 
PLS formation, since it is present in maximum amount in solutions de- 
rived from isolated paramyosin fibrils, and only in very small amount in 
“citrate” solutions which yield PLS most abundantly. Furthermore, PLS 
can be obtained from ultracentrifugal supernates of runs in which the 
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second component has been allowed to sediment out. The standard 
sedimentation constant (53) extrapolated to zero protein concentration at 
ionic strength 0.1 is about 4.3 Svedberg units. 

In electrophoretic analysis, the solutions exhibit a high degree of homo- 
geneity at ionic strength 0.1 or higher. However, the peak becomes in- 
creasingly asymmetrical and the pattern heterodisperse as the ionic 
strength is lowered. Diffusion runs under similar conditions (u & 0.1, 
pH = 3.5) give values for the diffusion constant of about 1.0 K 10 
Molecular weights computed from the sedimentation and diffusion data 


‘ 


are of the order of several hundred thousand, varying both with ionic 
strength and protein congentration. On this basis, utilizing the viscosity 
data, the calculated particle lengths lie in the range 1000-1500 A., the 
widths in the range of 15-30 A. Preliminary light scattering observations, 
employing both turbidity and disymmetry ratio determinations, confirm 
the presence of particles of similar molecular weight and length. Thus, it is 
clear that the major component of native paramyosin exists in acid solu- 
tions in the form of very thin highly elongated particles. 

As the ionic strength of the acid solutions is reduced to very low values, a 
characteristic change takes place in the physical properties of the elongated 
particles. The boundaries in the ultracentrifuge, and in electrophoresis 
and diffusion, become increasingly asymmetrical, suggesting that the pro- 
tein is undergoing a dissociation under these conditions. The asymmetry 
is always such that the peak spreads more toward the protein side of the 
boundary. The diffusion constant (uncorrected for asymmetry) increases 
by a factor of about 10. The electrophoretic pattern shows considerable 
inhomogeneity, of a type consistent with dissociation of a complex into a 
multiplicity of components having a wide spectrum of mobility characteris- 
tics. This suggests that the dissociation involves the breakdown of a com- 


plex rather than a simple depolymerization. The intrinsic viscosity does 


not undergo a drastic change as a result of the dissociation but rather 
increases by a relatively small factor. 

The physical chemical properties of the solutions seem to point to the 
following tentative conclusions: 

(a) The apparently homogeneous major component found in media of 
relatively high ionic strength (order of 0.1), and which in itself seems ca- 
pable of giving rise to PLS, appears to represent a complex of a highly 
asymmetric component and a second uncharacterized component. 

(b) PLS fibrils arise from these elongated complex units by a process of 
lateral aggregation as the ionic strength is increased. 

A more detailed physical chemical analysis 1s proceeding and the results 


will be reported in due course. 


* The data presented here are derived from a Ph.D. thesis in the Biology Department, 
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Massachusetts Institute of Technology, June, 1952. It is a pleasure to acknowledge the 
guidance and encouragement of Prof. F. O. Schmitt in this work 
t The author is indebted to the Commonwealth Scientific and Industrial Research 
Organization, Melbourne, Australia, for a Senior Studentship during the tenure of which 
much of the work described here was carried out. The investigation was supported in 
part by a grant from the Trustees under the wills of Charles A. King and Marjorie King 
' Hall, C. E., Jakus, M. A., and Schmitt, F. O., J. App. Phys., 16, 459 (1945) 


2 Bear, R. S., J. Am. Chem. Soc., 66, 2043 (1944). Cannan, Cecily, M. M., Ph.D 


Thesis, MIT, June, 1950 

3 However, the small angle diffraction pattern of intact clam muscle shows, in addi 
tion to the paramyosin reflections, a set of reflections corresponding to the so-called 
“myosin” net. The latter reflections are relatively strong in the red and weak in the 
white portions of the muscle, an observation in accord with the fact that the white por- 
tions contain a larger proportion of paramyosin. Furthermore, it seems likely that the 
treatment used in preparation of purified suspensions of paramyosin fibrils (maceration 
in dilute salt solution and repeated centrifugation) results in the elimination of the ma- 
terial giving rise to the ‘‘myosin”’ net reflections (presumably actoimyosin } 

* Highbeiger, J. H., Gross, J., and Schmitt, F. O., Proc. Nari. Acap. Sct., 37, 286, 
1951 


THE INFLUENCE OF DIPOLE MOMENT FLUCTUATIONS 
ON THE DIELECTRIC INCREMENT OF 
PROTEINS IN  SOLUTION* 


By Joun G. Kirkwoop AND JouN B. SHUMAKERT 
STERLING CHEMISTRY LABORATORY, YALE UNIVERSITY, NEW HAVEN, CONNECTICUT 
Communicated July 23, 1952 


Proteins in solution are known to produce! remarkably large increments 
in the dielectric constant of water. These increments have been con 
ventionally attributed to large permanent molecular electric dipole mo 
ments, which are partially oriented by an external electric field. An 
alternative mechanism of dielectric polarization of macromolecules con 
taining a large number of loosely bound ions, for example, protons, is 
provided by the migration of the small ions within the molecule under the 
action of the external field. For proteins, containing a number of neutral 
and negatively charged basic sites, such as = NH» and — COO™, to which 
protons are bound, this mechanism might be especially significant. Except 
in highly acid solutions, the number of basic sites generally exceeds the 
average number of protons bound to the molecule, so that there exist 
many possible configurations of the protons, differing little in free energy. 
Fluctuations in the configuration of the protons should be associated with 
a non-vanishing mean square electric dipole moment, even if the mean 


permanent moment were zero. In the present article we shall investigate 
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the contribution to the dielectric polarization of protein molecules arising 
from fluctuations in charge and configuration of mobile protons in the 
molecule. We shall show that the dielectric increments of many proteins 
can be entirely accounted for by this mechanism of polarization without 
postulating the existence of permanent electric moments. Fluctuations 
in the charge and configuration of mobile protons also make an important 
contribution to the force bet veen protein molecules. A theory of inter- 
molecular forces of this new type will be presented in a later article. 

For liquid solutions of high dielectric constant, the theories of Onsager* 
and Kirkwood? lead to the following relation between the dielectric con- 
stant D and the molar polarizations P, of the several components, 


9 cn 
D=1+5)P,; 


lorN 


se 


ce 


P; la, + (u,” fav. 3kT’} (1) 
where a, is the optical polarizability and (u,;?),y. is the mean square of 
the dipole moment of a molecule of type 7. Correlations in the orientation 
of neighboring molecules, treated in the theory of Kirkwood, will Le 
neglected here, an approximation adequate for protein molecules in 
aqueous solution. For the mean square dipole moment of a protein 
molecule, averaged over-all configurations of its mobile protons, we shall 
write, 
(uy? av x (u av + Ap’? 


Ap? = ((u = (u av lav. (2) 


where (u)ay. is the mean permanent dipole moment and Ay? is the dipole 
moment fluctuation associated with fluctuations in charge and configura- 
tion of the mobile protons. 

We consider a protein molecule cofitaining v basic groups of intrinsic 
charges, e;, situated at points R, relative to the molecule center of mass. 
We define a proton occupation variable x, for each basic group to be unity 
when that group is occupied by a proton and zero otherwise, The electric 
moment pw of the molecule, its mean value (u)sy, and fluctuation Ay? may 
then be expressed in the form, 


(u AV. >} (é; + e (xX; av Ry 
’ 1 


v 


Ap’ S e° Hi [ (xy, av. 


t,k= 1 
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where e is the protonic charge, and the averages are to be taken over all 
distributions, x;, -- - x, of the protons among the basic sites of the molecule. 
It is to be remarked that nesther (u),,. nor Aw? are indepencent of the origin 
to which they are referred, even at the isoionic point, because of fluctua- 
tions in net charge. They should properly be referred to the center of 
drag of the molecule, which is here assumed identical with the center of 
mass, which is true if non-linear hydrodynamic terms are neglected. 


1 
+ . BIA 
yx 


(4) 


where W’,(x;.. .x,) is the local free energy of proton configuration, x. ..x,, 
and A, is the total configurational free energy. The pertinent mean values 
(X; av. ANd (xix ay, OVer-all proton configurations may be calculated by 
methods previously developed by Kirkwood! in the theory of acid-base 


equilibrium of ampholytes. 
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where + is the radius of the protein molecule, assumed spherical, Ay is 
the intrinsic dissociation constant of the acid conjugate to the basis group 
i, [Z7*] ic the hydrogen ion activity, and « is the ionic strength parameter 
of the Debye-Hiickel theory. If electrostatic interaction between the 


protons is neglected, we obtain the simple result, 
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l 
1+ Ky/|H*} 
1#k 
l 
2+ K,/[Ht| + (H*)/K, 
In this approximation, the dipole moment fluctuation becomes 


pt ; 
Ay? =e 2 2+ K,/(H+*] + (H+//K, Me 


1 | = 
If there are v, equivalent basic groups of type @ in the molecule we may 
write, 
_ Vala” 
Au* = e* : : (S 
‘ sy K,/{H*+|] + [H+)/K, 


where KX,’ is the mean square distance of groups of type a from the molec 
ular center of mass. A further simplification is possible, if we assume 
that all basic groups are uniformly distributed on the surface of the mole- 


cule, considered as a prolate ellipsoid of revolution. 


Vv 


a 


Au? = e*ftby? D 
Ht ofibe 2y OK, [*) + (VK, 


+ 2)V o* 1+ o*(o* + 4) sec. ' a 


a’V ao 1+ o'sec.! a 


ach by (ab?) 


where o is the axial ratio of the ellipsoid and by is the radius of the equiv 
alent sphere. 

We present in table | calculations of dipole moment fluctuations, based 
on equation (9), for some representative proteins. It will be remarked 
that only those groups with pA values nearly equal to the pH of the solu 
tion contribute significantly to Ay’. At the pH values in question for 
the four proteins, @-lactoglobulin, ovalbumin, horse hemoglobin, and 
human serum albumin, only the carboxyl and histidyl groups make ap 
preciable contributions. For each protein, two values of Ay’ are pre 
sented, one for the molecule considered as a sphere, and one for the mole 
cule considered as an ellipsoid. The values are compared with (uy? beeen 
calculated from Oncley’s! measurements of the dielectric increments of 
the several proteims. The experimental values are subject to considerable 
uncertainty and will doubtless undergo revision on the basis of more pre 
cise measurements. Nevertheless, it will be observed that the calculated 














cr nana 





VoL. 38, 1952 CHEMISTRY: KIRKWOOD AND SHUMAKER 859 


values of Aw are of the same magnitude as the experimental values of 
(2), /*. Only in the case of 8-lactoglobulin does it seem to be necessary 
to postulate the existence of a permanent dipole moment in order to account 
for the observed value of (u?),,."’. The dielectric constant increments 
of ovalbumin, hemoglobin and serum albumin are all adequately accounted 


rABLE 1 


Dipole MOMENT FLUCTUATIONS WitHouTr ELECTROSTATIC INTERACTION 


HORSE HUMAN 

B-LACTO OVAL HEMO SERUM 

PROTEIN GLOBULIN BUMIN GLOBIN ALBUMIN 

Molecular weight 37,000 16,000 68 ,O000 69,000 

Dissociating groups per molecule“ 

Arginine (pK = 12.5) 6 15 14 25 
Histidine (pk = 6.0) i 7 36 16 
Lysine (pK = 10.5) 20 20 38 58 
Tyrosine (pK = 10.1) 8 9 11 18 
Free carboxyl (pK = 4.0) 59 51 53 93 
pH 5.6 1.8 6.4 5.0 
by (A.) 24 26 30 30 
a/b 5 1 5 5 
(u?)ay. ? Debye units (experimental ) 700 260 480) 380° 
Au (sphere ) 170 320 390 130 
Ap (ellipsoid ) 270 140 620 680 


“CG. R. Tristram, Advances in Protein Chem., 5, 84 (1949) 
’ This experimental dipole moment is for horse serum albumin. Professor J. L. 
Oncley has informed us that recent measurements of the dielectric increment of serum 
. > : » i 
albumin lead to a value of 700 Debye units for (u?)av 
TABLE 2 


DreoLE MOMENT FLUCTUATIONS WITH ELECTROSTATIC INTERACTION 


B-LACTO OVAI HEMO SERUM 
PROTEIN GLOBULIN BUMIN GLOBIN ALBUMIN 
pH 5.6 18 6.4 5.0 
Au (sphere ) 360 330 300 320 
Ap (ellipsoid ) 560 160 180) 510 
rABLE 3 


Errecr OF PH ON THE DIPOLE MOMENT FLUCTUATION OF OVALBUMIN 

pH 2.0 1.0 18 6.0 8.0 10.0 

Am. : 149 390 330 224 163 5 
for on a semiquantitative basis, by the calculated values of Ay? arising 
from fluctuations in charge and configuration of the mobile protons. 

The effect of electrostatic interaction between the protons has been 
calculated numerically at zero ionic strength with the use of equation (5) 
in its complete form, for the four proteins. The results are presented in 
table 2. It will be remarked that Ay is closer to the experimental value 


of (u®),..* for B-lactoglobulin than without interaction. Although the 
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FIGURE 1. 

The effect of pH on the dipole moment fluctuation of ovalbumin without electrostatie 








interaction. 
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FIGURE 2. 


The effect of pH on the dipole moment fluctuation of ovalbumin with electrostatic 
interaction. 
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influence of interaction is by no means negligible, it does not affect the 
magnitude of the dipole moment fluctuation. Electrostatic interaction 
causes configurations of low net charge to be favored and is most important 
at values of pH departing from the isoionic point. 

From equations (3) and (4), we obtain the following expression for the 
mean dipole moment, with the neglect of electrustatic interaction between 


the protons, 


eR, 


= 


= _ 
” r Feast l 13 K, (i+ | 


+ 


w= > eR, (10) 
i ] 

where w is the permanent dipole moment of the fixed charges of the mole- 
cule, exclusive of that of the protons. Unless each type of basie group 
is distributed on the surface of the molecule with a center of symmetry, 
uw does not vanish and (w),y cannot vanish at all pH. Indeed at high 
values of pH, (w)av. approaches wo. If wo is opposite in sign to the second 
terms of equation (10), as would be expected to be the case if all basic 
groups are either neutral or negatively charged, the absolute magnitude 
of (u)ay. would be a monotonically increasing function of pH. The dipole 
moment fluctuation, Au’, on the other hand, possesses maxima at pH values 
coincident with the pK of the conjugate acids and minima at intermediate 
values, as illustrated by figure 1, for ovalbumin. Away from the isoionic 
point, the effect of electrostatic interaction suppresses the maxima and 
minima in Ay’, as illustrated for ovalbumin in table 3 and figure 2. It 
will be observed that Ay® is a monotonically decreasing function of pH, 
diminishing from a high value at pH 2 to a value effectively zero at pH 10. 
The experimental data relating to the effect of pH on the dielectric incre- 
ment of Shutt® and of Oncley, Bateman, Pecher and Melin' are too frag- 
mentary and contradictory to justify an attempt at quantitative inter 
pretation by means of the present theory. Since the effect of pH provides 
an important criterion for determining the relative importance of the con- 
tributions of (u)ay.? and Au’ to the dielectric increment, additional experi 
mental data would be highly desirable. 

The theory which has been presented here relates only to the dielectric 
merements of proteins under the action of static electric fields. A theory 
of dielectric dispersion and absorption by the proposed mechanism of 
fluctuations in proton configuration remains to be developed. The re- 
laxation time spectrum of an ellipsoidal molecule will be determined not 
only by external rotatory diffusion but also by the diffusion of the mobile 
protons on the surface of the molecule. This means that the structural 
interpretation of dielectric dispersion will be considerably more com- 
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plicated than for a molecule possessing a permanent dipole moment with 
out fluctuation. Qualitative considerations, however, lead to the con- 
clusion that internal relaxation times of the order of magnitude of those 
observed are not inconsistent with the proposed mechanism. For example, 
for a sphere of radius }, on the surface of which protons possess a trans 
latory diffusion constant D, a relaxation time of the order of magnitude 

b’?/D should control the distribution of protons. With 6 equal to 
30 A., the radius of a sphere of volume equal to that of a serum albumin 
molecule and with D equal to 10° em.*/sec., a relaxation time of the 
order of 10° sec. is obtained, which agrees in order of magnitude with 
those observed for representative proteins. 

According to the present theory, the structural significance of the 
dielectric increments of proteins is subject to rather drastic reinterpreta- 
tion. While many proteins doubtless possess permanent electric dipole 
moments of considerable magnitude, they cannot be determined un- 
ambiguously from measurements of the dielectric increment, since the 
dipole moment fluctuation, arising from the mobile proton distribution 
has been shown to make a contribution to the increments, which is equiv 


alent to that of a permanent moment of large magnitude. The influence 


of fluctuations in the charge and configuration of mobile protons on the 
forces between protein molecules, which is of greater significance than the 
effect on the dielectric increment, will be treated 1n a later article. 


* This investigation was in part supported by the office of Naval Research 

t Part of a dissertation presented for the degree of Doctor of Philosophy in Yale 
University 

'Oncley, J. L., Chapter 22, Cohn and Edsall, Proteins, Amino Acids and Peptides, 
Reinhold, New York (1945) 

2 Onsager, L., J. Am. Chem. Soc., 58, 1486 (1936) 

* Kirkwood, J. G., J. Chem. Phys., 7, 911 (1989). 

‘ Kirkwood, J. G., Chapter 12, Cohn and Edsall, Proteins, Amino Acids and Peptides, 
Reinhold, New York (1943). 

®Shutt, W. J., 7rans. Faraday Soc., 30, 893 (1934) 
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FORCES BETWEEN PROTEIN MOLECULES IN| SOLUTION 
ARISING FROM FLUCTUATIONS IN PROTON CHARGE 
AND CONFIGURATION* 


By JoHN G. KIRKWOOD AND JOHN B. SHUMAKER 


STERLING CHEMISTRY LABORATORY, YALE UNIVERSITY, New HAVEN, CONNECTICU1 


Communicated July 23, 1952 


The forces between protein molecules in solution are in large part elec- 
trostatic in origin. Evidence for the dominant role of electrostatic forces 
is provided by the marked effect of ionic strength on the thermodynamic 
interaction of proteins. Thus, the commonly observed reduction in 
interaction produced by the screening action of the statistical space 
charge of the 1ons of an electrolyte would only be effective on that part 
of the force between protein molecules, which is electrostatic in origin. 
At values of pH departing from the isoionic points of the molecules in 
question, a simple Coulomb force, determined by their net electric charges, 
is predominant. This force is non-specific and is sensitive to molecular 
structure only in so far as structure influences the net charges. When one 
or both of a pair of protein molecules are at their isoionic points, structure 
sensitive electrostatic forces come into play. These forces have been 
customarily attributed to fixed constellations of electric charge, which 
impart to the molecules permanent electric multipole moments charac- 
teristic of their structure. 

We shall investigate here another type of electrostatic interaction 
between protein molecules, which arises from fluctuations in charge and 
charge distribution associated with fluctuations in number and con- 
figuration of the protons bound to the molecules. Proteins, considered as 


ampholytes, contain a large number of neutral and negatively charged 


basic groups, such as — NH» and ~~ COO, to which protons are attached. 
Except in highly acid solutions, the number of basic sites generally exceeds 
the average number of protons bound to the molecule, so that there exist 
many possible configurations of the protons, differing little in free energy 
among which fluctuations may occur. Fluctuations in the number and 
configuration of the mobile protons impart to the molecules fluctuating 
charges and fluctuating electric multipole moments. We have demon 
strated in a previous investigation! that the dipole moment fluctuations, 
arising from configurational fluctuations of the protons, are sufficient to 
account for the dielectric constant increments of many proteins without 
postulating the existence of permanent dipole moments. Linderstrém- 
Lang? has shown that fluctuations in net charge contribute to interaction 
between proteins and the small ions of an electrolyte. In the present 
investigation, we shall demonstrate that the interaction between the 
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fluctuating charges and multipole moments of two protein molecules 
makes a significant contribution to the intermolecular force, by a relatively 
simple mechanism. The fluctuating electric field of each molecule alters 
the distribution of fluctuations in the charge and constellation of the 
mobile protons of the other in such a manner as to produce at the isoionic 
point a long range attractive force between them with a potential diminish- 
ing asymptotically as 1/R’, in the absence of screening by the statistical 
space charge of an electrolytic environment. With screening, the range 


«R 


is reduced by an exponential factor, e “", depending upon ionic strength 
through the Debye-Hiickel parameter, x. 

We consider two protein molecules in aqueous solution, the centers of 
mass of which are separated by a distance R. We suppose that molecule 
| contains », basic groups of intrinsic charge e,"’) and that molecule 2 
contains v, basic groups of intrinsic charge e,. We define proton oc- 
cupation variables x," and x,“ to be unity when the respective basic 
groups are occupied by one proton and zero otherwise. For specified 
proton configurations, x") x, and x, ....x,,", the electric charges 
of the several groups, q,") and gq are e? + ex; and e + ex,, 
where e is the protonic charge. The mutual electrostatic energy V of the 
two molecules, in specified orientations and proton configurations is given 
by, 


~ |" 
. = a 


qi‘? - e,? t ex," 


qe? -_ e,°?) + ex, (2) (1) 
where D is the dielectric constant of the solvent and R,,“'” is the scalar 
distance between basic group 7 of molecule | and basic group k of molecule 
2. Electrolytic screening is not included in equation (1), but will be 
taken into account for special models later. For generality, the sum of 
equation (1) may be regarded as extending over other charged groups of 
the molecules, as well as basic groups, with the convention the proton 
occupation variable is always zero for a non-basic group. 

The potential, W(R), of average force between the two molecules is 
related to the potential |’ of the force in fixed orientation and proton 
configuration in the following manner, 


BWR 
oe BW(R) 


W(R) = eve = P Pay) ee OCB) 





Vor. 38, 1952 CHEMISTRY: KIRKWOOD AND SHUMAKER 865 


where the average is to be taken over all orientations and over all proton 
configurations of the two molecules, with uncorrelated distribution func- 
tions. For simplicity in the present argument, we shall suppose that at 
distances not permitting overlap of the excluded volumes of the two mole- 
cules, their only interaction is electrostatic, with the potential, equation 
(1), and we shall retain only the first two terms in the expansion, equation 
(2) of W(R) in powers of 8. With these simplifications, we may express 
W(R) as follows, 
(gq), <q), no 
wim) = =. = * — 2» 1 » 1 Oa Age aw S 
R? 
Ff xis = NR? Ry? gy. (3) 


where (gay. and (g™),y. are the average total charges of the two 
molecules, equal to zero at their respective isoionic points, and the func- 
tions fun, averages over all orientations of both molecules, approach 
unity asymptotically as the distance R between their centers of mass 
becomes large relative to the linear dimensions of both molecules. We 
now write, 

(gi? )av. + Age” 

(qe av. + Aq” 

efx? — (x, ay. | 

= e(x,‘” (x, 


(4) 


in order to separate the contributions of the mean and fluctuating charge 
distributions to W(R) in the following manner, 


WR) = Fee. 4 yyencR) + WOR) + WD(R) + WO(R) 
DR 


B a 
W)(R) = (¢ i” av (¢ Da (qx? av .\4 AS» ie s 
oDIR? 2, f] 11 1k 1 fii 


WO (R) 7 ‘ “ 9 (Agi? Agi av. (ge) av (qa) av S ints 


W°YO(R) = — (qa dav. (Ga? dav (Age? Age av. S sets 


VOD(R) a (Agi Agi) av (Agu Aq, ae S tuts 
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By expansion of the functions /,,, in powers of (1/R), each of the terms 
W'*”(R) may be expressed as a sum of terms corresponding to the inter- 
action of electric multipoles of various orders. ‘The term W(R) repre- 
sents the interaction of the permanent electric multipoles of the molecules, 
the terms W°"(R) and W"(R) the interaction of the permanent multi- 
poles of the one molecule with the fluctuating multipoles of the other, and 
the term W°"'(R) the interaction of the fluctuating charges and fluctuating 
multipoles of both. 

We shall now limit our considerations to values of pH corresponding 
to the isoionic point of one of the proteins. Under this condition, W(R) 
becomes asymptotically equal to the fluctuation term, W'!'(R), at large 
intermolecular distances, as each of the functions fi,;, approaches unity. 
Among the several contributions to W(R), we shall therefore undertake a 
detailed analysis only of this term. With the use of equations (4) and 


(5), it may be expressed in the form, 
WOO(R) = - 
pe Rk? om xk 
Ay CAL hay If iets (6) 
The pertinent mean values of the products of the proton oo vari- 


ables are to be taken over distributions, x,"' ye ND: opel Ch 
of the protons among the basic sites of the two molecules. For wads ip 


we have, 


x1 veh 


where WW," is the local free energy of proton configuration x") 
and A,‘" is the total configurational free energy of molecule |. The 


configurational mean values for molecule 2 may be calculated from exactly 
similar relations. Methods developed by Kirkwood*® in the theory of 
acid-base equilibrium of ampholytes may be employed for the explicit 
calculation of the mean values of the occupation variable products. 


0 log G"' 
0 log XA," 


(Xj dav 


A PA,O OG" 
G) Ad,MOA,4 
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B(x, 


2b,;DkRT l t Kh, (S) 


where },; is the radius o: molecule | considered as a sphere, A;‘!’ is the 
intrinsic dissociation constant of the acid conjugate to the basic group 2, 
|H7+| is the hydrogen ion activity, and « is the ionic strength parameter 
of the Debye-Hiickel theory. The averages for molecule 2 may be cal- 
culated by means of similar relations. If intramolecular electrostatic 
interaction between the protons of each molecule is neglected, the following 
result is obtained, 


ii 
2+ Ky/(H*)] + (Ht)/KY" 


e 


Oks 


2+ K,/(H+*] + (H*)/Ki@ 


(9) 
where 6,,, the Kronecker delta, is equal to unity of 7 / and zero other- 
wise. In this approximation the potential H''(R) becomes 


Be ‘ 


W'Y(R) 


ys > [2 + K,/[H*] + 
I 


. open = 


(I+)/K,]— [2 + Ky /[H*] + [H+)/K, 
| F l (u, 4 =A 
Og : 
hu ju, | l (u; Uy) 
(1 4 u,)(1 + ty — Uy) | 
u, log 
(1 + uy; uz)(l — uu; U,) 
; Uy)(1 + ty — Uy) 
u, log ( 
(1 - — u,)(1 — uy — Ux) 
u; 6 " D/R 


1 


where 6," and b, are the distances of the basic groups 7 and & from the 
centers of mass of the respective molecules | and 2, and the symbol for 
the function f,,;, has been abbreviated to fy. The factor fy, does not de- 
part greatly from unity. For a pair of identical spherical molecules with 
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basic groups located on their surfaces, /,, diminishes from 1.39 on contact 
to its asymptotic value, unity, at large separations. The asymptotic 
form of W°"(R) is obtained by approximating all /,, by unity. 
(Aq")av (Aq? ar 
WO (R) el A oe 
2D°R°k1 


(Ag? )a, 


(Ag? yw. =e DO (H+) + (H*V/K,@]-! (11) 


where (Aqg"! 7 vo: eAtiGl (Ag’?” wy. are the total charge fluctuations of the two 
molecules. 

The potential W""(R) of equation (11) leads to a long-range attractive 
force between the two molecules, diminishing as the inverse cube of the 
distance. A force of this range fails to yield convergent expressions for 
the thermodynamic functions. Properly, a Debye-Htickel factor should 
be included to provide for screening of the fluctuating electrostatic inter- 
action by the statistical space charge arising from counterions to the pro- 
teins and from the ions of other electrolytes present in the solution. If 
the molecules are spherical, the asymptotic potential, equation (11), is 
modified by screening as follows, 


2(R dir) 


. (Aq! aN (Ag’?”’ ke 
WOPCR) ~ 
2D?R°k1 (1 + Kay)? 
Ko" { Ki" 
inN F(Ag vay Cy : (Ag’?”)ey 4] 
ioopeT | vA VV 
lr Ne? - 
LOOODRT 


(12) 


where dj is the sum of the radii of the two molecules, x’ is proportional 
to the ionic strength of other electrolytes present in the solution and the 
second term represents the contribution of the fluctuating charges of the 
proteins to «*, their concentrations C; and C, being expressed in grams 
per 100 ml. of solution. Here A/; and Mp», are the molecular weights of 
the proteins and N is Avogadro's number. If either of the proteins is 
not at its isoionic point, (Ag®)ay. should be replaced by the mean square 
<q’> ay, to its total charge. When either of the proteins is isoionic, 
WD(R), of equation (12), is the dominant asymptotic term in the total 
potential of average force, IV(R), at large intermolecular distances, since 
then the first term of equation (3), appropriately modified with a Debye- 
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Hiickel screening factor vanishes. A more exact calculation of W"(R), 
based on equations (6), (7) and (8), taking into account the influence of 
intramolecular electrostatic interaction between the protons on the mean 
values of the proton occupation variable products, may be carried out 
numerically. When both proteins are isoionic, the correction for proton 
interaction is not large. 

According to the theory of Kirkwood and Shumaker,' W“"(R) may 
also be expressed in terms of the dipole moment fluctuations Ay,’ and 
Ap,’ of the two molecules, through their relationship to the charge fluctua 
tions. The expression is the following, 

‘ . ®(R ly 
ORE). see: (13) 
2D°b2 be kT R*(1 + Kay2)? 
where }, and /, are the radii of the two spherical protein molecules. Equa- 
tion (13) is convenient for estimating the order of magnitude of W°"(R) 
from experimental values of the dielectric increments of the proteins. 

We shall now investigate the influence of the interaction of protein 
molecules through charge fluctuations on the chemical potential of a single 
protein at its isoionic point. According to the theory of Kirkwood and 
Buff,4 we may write at low concentrations, for a protein component | in 
a solvent 0, at constant xo of other electrolytes present in the solution, 


LOOM, (Ou. ‘ . 
on de = Gy ~ Gi 
Nk7 Oc, / 1 DX, 


Gy lor / R*{gy(R) — |] dR 


Gio i f R? g( R) . 11dR 


My kT log Cy + wie + wi(7, pP) 


mi(7’, p) lim [u; — kT log Ci} (14) 


c->0 


where g1,(R) and gy(R) are the radial distribution functions of a pair of 
protein molecules and a protein and solvent molecule, respectively, which 
are related to the potentials of Average force Wi,(R) and Wyw(R) in the 
following manner, 

g1,(R) e BWiiR) 

gu(R) = e~ OW) (15) 


In each of the integrals of equation (14) there is a co-volume term arising 
from the region of overlap of the excluded volumes of the molecular pair. 
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If the dimensions of a solvent molecule are ignored in comparison with 
those of the protein, the co-volume contribution to Gy is one-eighth that 
of Gy. If, further, the solvent is assumed to play only the role of a di- 
electric continuum exterior to the excluded volume of the protein, and 
W,(R) is approximated by (12) exterior to this volume, we obtain with 
the neglect of powers of order 6’, 

LOOM, (™) on wN(Aqg’)ay 

kT \de/t. 0,5 ~~ (DRT) *«(1 + wa)? 
B, = 7xNa*/12 (16) 

where a is the molecular diameter and (Ag‘")ay. is given by equation 
(11). If the contribution of the protein to the ionic strength is negligible 


in comparison to that of other electrolytes present in the solutions, we 
obtain by integration of equation (16) the following result for the excess 


chemical potential, u,, per molecule. 


mi /kT = 2BC, 


s=( 3 ip mN(Ag)ev. (17) 
~ \00M SU? (DRT) %xo(1 + Koa) 28 P 


The coefficient of C,; in the expansion of the osmotic coefficient in powers 
of C; is equal to B by conventional thermodynamic relationships. In the 


absence of other electrolytes with the ionic strength determined entirely 
by the fluctuating charge of the protein, integration of equation (16) 
yields the Debye-Hiickel limiting law, 
ecLmT (Aq )ay 9 
My /k7 - —_ DRT Ki -+ OK, ) 
tr( Ag’ av. 
lOO.M,DRT 


characteristic of electrolytes. 

At high ionic strengths, the co-volume contribution, By, to the inter- 
action cveflicient, B, of equation (17) is the dominant one leading to a 
positive value of the excess chemical potential. Only at quite low ionic 
strengths does the second term arising from proton fluctuations become 
dominant, leading to a negative excess chemical potential. For serum 
albumin in water at 300 K at its isoionic point, the coefficient B is estimated 


1.2\ / 
as (12 Se )/ 1000 (19) 
VIy/ 


/ 


to be, 


where I’ is the ionic strength of other electrolyte. At ionic strength of 
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the order .f 10~*, the fluctuation term is approximately ten times as large 
as the co-volume term and gives rise to substantial negative values of the 
excess chemical potential at protein concentrations of the order of several 
per cent. 

The present analysis of the attractive force between proteins arising 
from fluctuations in charge and configuration of mobile protons is neces- 
sarily schematic because of lack of knowledge of the details of protein 
structure. Although it is a force of long range at low ionic strength, it 
appears to exhibit specificity only through the influence of structure on 
the fluctuations. It is clear, however, that highly specific interactions 
might well arise from the mechanism of interaction, which has been de- 
scribed. In favorable orientations, steric matching of a constellation of 
basic groups on one molecule with a complementary constellation on the 
other could conceivably produce a redistribution of protons leading to a 
strong and specific attraction depending upon the local structural details 
of the complementary constellations. As an extreme example it is possible 
to imagine proton fluctuations to be frozen in such a manner as to produce 
an intermolecular zwitterion, with matching areas of positive and negative 
charge on complementary areas of the two molecules. Such considerations 
relating to specificity of the fluctuation force must necessarily remain 
speculative until detailed knowledge of structure is available. 

In conclusion, the authors wish to acknowledge their indebtedness to 


Julian M. Sturtevant for interesting discussions during the course of this 


investigation. 
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AN ASSOCIATION OF HEREDITARY EYE DEFECTS WITH 
WHITE SPOTTING 
By L. C. DUNN AND JAN Monr* 


CoLUMBIA UNIVERSITY 


Communicated August 20, 1952 


In a paper in these PROCEEDINGS! an asymmetrical effect on eye color 
in the house mouse was ascribed to a new mutant gene, ruby eye (ru). It 
was noted, however, that almost all mice in which the two eyes differed 
in color, referred to as heterochromia .ridis, showed piebald spotting. This 
connection between eye color variation and spotting has now been studied 
further with the result that the eye color modification has been shown to 
be due to specific defects of the iris, which in turn appear to be effects of a 
mutant gene, s, for piebald, which has long been known in laboratory mice. 
A strain of ruby-eyed mice without spotting has been observed for over 
twenty generations and no case of heterochromia has been seen. 

Evidence for the association of iris defects and spotting will be given be- 
low, together with a preliminary description of the iris abnormalities de- 
rived from gross examinations with an ophthalmoscope. All animals in 
these experiments were thus observed at the time of weaning (3-4 weeks of 
age) and the diagnosis was in many cases confirmed by later examinations. 
A careful microscopic study of the defects remains to be made. 

Anisocoria in Piebald Mice..-Many of the mice in our inbred strain 
An (= aniridia: black-and-tan silver ruby piebald a‘a' BB ruru sist ss) 
show noticeable differences from the normal ruby eye color. One or both 
eyes may be light ruby or pale pink. When these variant eyes are ex- 
amined with an ophthalmoscope the commonest finding is an enlarged 
pupil (anisocoria) varying all the way from slight enlargement to what 


appears to be complete absence of the iris (aniridia).. Enlargement of the 


pupil is generally accompanied by reduction of the amount of pigment in 
the iris with consequent approach to the pink color of unpigmented eyes. 
A smaller proportion of lighter colored eyes have pupils of normal or nearly 
normal size but lack a part or all of the dark iris pigment. Occasionally a 
segment of iris is unpigmented, Out of 266 eyes of this strain examined in 
mice of the sixth and seventh inbred (brother-sister) generations, 120 had 
the iris reduced (98) or missing (22); 4S without gross iris defects had pig- 
ment reduced (34) or missing (14); while the remaining 9S had no gross 
defects. Of the 133 animals examined, both irides were normal in only 28, 
so that in this strain the iris defect expresses itself in about SO per cent of 
the individuals and in some 63 per cent of the eyes. Thus there is a high de- 
gree of asymmetry in gross expression. The right and left eyes have about 
the same probability of being abnormal as judged by ophthalmoscopic ex- 
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amination. The proportion of eyes with histological abnormalities is 
unknown. The differences in iris character between two sides of the same 
individual and between individuals of this inbred strain are apparently 
non-genetic. ‘ 

In contrast with this strain is another of similar derivation. This strain 
(Anophthalmic), now in its 22nd brother-sister generation, has a genotype 
which is similar to that of the Aniridia strain except that it lacks the pie- 
bald gene (a‘a' BB ruru sisi). It contains a gene or genes which express 
themselves in some 65 per cent of the individuals, often asymmetrically, by 
causing a defect of the bulbus and/or cornea which varies from bilateral 
anophthalmia to unilateral staphyloma. In none of the eyes of this strain 
with normal corneae have the iris defects of the Aniridia strain been seen. 

Tests of Inheritance of Anisocoria.—The most direct test of the hypothesis 
that anisocoria is an effect of the piebald gene is to cross the above two 
strains which differ primarily at the piebald locus. From this cross 26 F; 
animals were examined. All were non-piebald with ruby eyes and normal 
irides. 

Fourteen of these F;’s were backcrossed to the spotted ruby (An) strain 
and produced the results shown below: 


NORMAL IRIDES ABNORMAL IRIDES TOTAL 
Non-piebald 116 } 120 
Piebald 72 13 1L5 


Of the non-piebald offspring only 3.3 per cent had any abnormality of 
the iris, and of these one had an oval pupil on one side, and one had a 
slightly enlarged pupil on one side. Of the other two, one had no detect- 
able iris on one side; the other had an unquestionably enlarged pupil on 
one side. These last two cases are exceptions to the rule that anisocoria 
appears only in piebald animals. ‘They may of course have been homozy- 
gous for the piebald gene ss which for some reason did not express itself in 
spotting, but there is no experimental test of this possibility. Two cases 
of what then was called heterochromia iridis were noted among 48 non- 
piebald ruby eyed mice.' Although these were not examined with the 
ophthalmoscope, they were probably exceptions of the sort noted above. 

The significant feature of the results was the appearance of iris abnor- 
malities in some 37 per cent of the piebald animals. This is not far from 
the frequency of iris abnormalities in the An piebald strain at the time 
(F;) when the cross was made (15 out of 52 animals examined showed 
gross abnormalities of the iris). 

The data show that piebald spotting and iris defects are closely asso- 
ciated and do not exclude the hypothesis that they are due to the same 
genotypic cause. 

Anitsocoria in Other Spotted Strains.-lf the iris defects are due to the 
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action of spotting genes they should be found in other spotted stocks un- 
related to the above. Among 39 piebald mice from an inbred black Brachy 
line in this laboratory known to be genetically ss, 5 cases of small or absent 
iris were found. Among 63 piebald animals af Line 3 (Black agouti) in 
this laboratory, eleven had gross defects of the iris. One case was found 
among 95 animals with spotted face in Dr. Schneider’s colony at Rocke- 
feller Institute for Medical Research and one case among 15 animals with 
the “splotch’ mutant condition in strain C57 black. The last two strains 
are not known to carry the piebald gene. These 1S cases were all found 
in black-eyed animals so that although the ruby-eye phenotype makes 
the iris abnormality easier to identify with the naked eye, ruby is not a 
necessary condition for its development. In all of the 102 black-eyed pie- 
bald animals from the Brachy line and Line 3 a deficiency of the chorioidal 
pigmentation could be recognized ophthalmoscopically. This confirms 
earlier findings of defective chorioidal pigmentation in piebald mice.** 
In ruby animals this deficiency is difficult to recognize ophthalmoscopi- 
cally, but histological examination indicates that it is a general feature of 
the piebald phenotype in these animals too. 

Test crosses between a piebald with abnormal iris from Line 3 and one 
from strain An produced 5 spotted animals with abnormal iris and 9 spot- 
ted without iris abnormality, showing a similar genetical basis for aniso- 
coria in those two unrelated strains. A test of one anisocoric animal from 
the Brachy strain by strain An was negative yielding only 15 animals with 
normal irides. Because of the small numbers involved, this test was not 
conclusive for the question at issue. 

Among the 176 self-colored animals from a variety of strains, examined 
with the ophthalmoscope, not a single case of atiisocoria was found. How- 
ever, 15 cases of abnormal iris were found among 100 albinos examined 
from various strains. Six of these cases, all from an inbred albino Fused 
strain of this laboratory, had coloboma iridis and were thus quite different 
from the anisocoria associated with piebald. The Fused strain was tested 
and found not to carry piebald. These and other cases of eye abnormalities 
found in several strains will provide good material for ophthalmological 
study. 

Finally, a new occurrence of Aniridia in several related animals from a 
mixed stock has been observed. In the first case, a brown piebald animal, 
both eyes were light pink. Wher tested by ruru SS this animal was shown 
to have the dominant allele of ruby (Ru), giving 26 non-ruby in F,, and 
matings among F, animals according to birth classification showed segre- 
gation, for uncolored irides (SO colored, 29 uncolored) for one eye color mu- 
tant only, viz., ruby eye. 

In F, and backeross progenies, ophthalmoscopic examination of irides at 
weaning revealed several cases of anisocoria including aniridia, but these 
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occurred only among spotted animals. Since the case was apparently 
similar to that analyzed above it was not studied further. It does show, 
however, that although anisocoria is easier to detect due to the pinkness of 
the eye, when it occurs in ruby-eyed animals, it may exert its full effect in 


spotted animals without the ruby mutation and thus simulate a mutation. 

Discussion. It is evident from the above that the type of related iris 
abnormalities first noted in strain An and shown there to be associated 
with the piebald condition are found also in other unrelated piebald strains 
of various genotypes and rarely in spotted strains not known to carry the 
gene S. 

This is probably not the first time that eye defects have been noted in 
piebald mice, although it is the first time that their genetic nature has been 
analyzed. Durham? in 1908 noted “ruby eyes” in three mice. In each 
case the choroid pigment was greatly reduced, and in each case the abnor- 
mal individual was chocolate (brown) pieba/d. No information on size 
of pupil was given. Later Gates* examined a population segregating for 
piebald, agouti, black-brown, and density and dilution, and found that all 
spotted animals and only spotted animals lacked chorioidal pigmentation. 
External defects, such as pinkness or heterochromia, were not noted. The 
present study indicates also that in piebald mice the chorioidal pigmenta- 
tion is generally defective. 

There is thus a general association between spotting and eye pigmenta 
tion which suggests the following working hypothesis: The genotype re- 
sponsible for piebald spotting (and perhaps for certain other spotting types) 
produces its effect by increasing the probability of death or failure in func 
tion of the pigment cells (melanophores). In the mouse these cells are de- 
rived from the neural crest (Rawles*) and reach their destinations by migra 
tion. The normal development of the iris depends upon actual or inductive 
contributions by such cells. If their probability of survival is lowered, the 
failure of pigment development in local areas of the skin will be accom 
panied by failure to invade the choroid, and an increased probability of 
failure of development of the stroma of the iris or of its pigment or of both. 
‘spotting gene’’ upon a struc- 


‘ 


In some such manner as this the effect of a 
tural component such as the iris may be envisaged. 

Since both amiridia and spotting occur as genetic variations in man, a 
study of their relations in human families may yield some evidence on the 
hypothesis suggested. 

* Visiting Fellow of the Rockefeller Foundation, 1949; now at Genetics Institute, 
The University, Oslo, Norway 

' Dunn, L. C., Proc. Nati. Acapb. Scr, 31, 343-346 (1945) 

? Durham, F. M., Reports to the Evolution Committee, Royal Society, 4, 41-57 
(1908) 

* Gates, W. H., Carnegie Institution of Washington, 337, 83-138 (1926). 

' Rawles, M. E., Physiol. Zool., 20, 248-266 (1947). 
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“BENT-TAIL,” A DOMINANT, SEX-LINKED MUTATION IN THE 
MOUSE* 


By E. D. GARBER 


NAVAL BIOLOGICAL LABORATORY, NAVAL SUPPLY CENTER, 
OAKLAND, CALIFORNIA 


Communicated by W. E. Castle, July 21, 1952 


Although sex-linkage in animals is no longer an unusual phenomenon, 
its apparent rarity in rodents is noteworthy, especially in the intensively 
studied house mouse and rat. Hauschka, ef al.,' have presented evidence 
for a sex-linked lethal altering the sex ratio in the house mouse; Falconer? 
recently reported a second sex-linked mutation in the house mouse. A 
third sex-linked mutation, Bent-tail (Bm), in the same species was indicated 
in a preliminary report.’ Additional data on the inheritance of Bent-tail 
are presented in this report. 


TABLE 1 


THE DISTRIBUTION OF FEMALES AND MALES ACCORDING TO THEIR EXPRESSION OF 
Bent-Tar (Bn) In VARIOUS TYPES OF MATINGS 


TAIL 
ONE-HALP 
TAIL NORMAL LENGTH ~ NORMAL 
CURLED, LENGTH, 
SINGLE SEVERAL SEVERAL SEVERAL 
SINGLE CLEAR CLEAR CLEAR CLEAR 
MATINGS NORMAL BEND BEND BENDS BENDS BENDS TOTAL 


Females 
bn bn 2 »& 
Bent-tail « 
Bn bn 9 X 
Bent-tail o 
Bn bn Q &X 
normal 


Males 
Bnbn 2 X 


normal 
Bnbn 2 & 
Bent-tail 


Origin and Description of the Mutation._-In the course of breeding mice 
for studies in experimental epidemiology, a single Bent-tail male occurred 
in a litter of 7 mice from a mating between a normal Namru strain‘ female 


hy 4 ° 
and a bald,’ fr’’, male. Subsequent litters were normal. 


The expression of Bent-tail in heterozygous females ranged from a single, 
almost imperceptible bend to the extreme cases in which the tail was curled 
and pressed against the body. The first 64 heterozygous females from 
matings 1 and 3 and the tirst 19 Bent-tail females (Bn bn and on Bn) 
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from matings between heterozygous females and Bent-tail males (mating 2) 
were scored according to an arbitrary classification of the expression of the 
mutation (table 1). The data from matings | and 2 are more reliable as 
regards expressivity than those from mating 3. All females in litters from 
mating | were heterozygous (Bn bn) and the females in litters from mating 
2 were either heterozygous or homozygous Bent-tail (Bn bn or Bn Bn); 
litters from mating 8, however, contained both normal and heterozygous 
Bent-tail females. The first 23 Bent-tail males from matings between 
heterozygous females and either normal or Bent-tail males were scored in a 
similar manner (table 1). The expression of Bent-tail in heterozygous 
females was more variable than in Bent-tail males. Penetrance in females 
was very high, at least 95°); penetrance in males was 100°. Tail length 
was normal in heterozygous females but approximately half normal length 
in Bent-tail males and in homozygous females, 


TABLE 2 
A SUMMARY OF RESULTS OF MATINGS INVOLVING BENT-TAIL (Bn) 
PROGENY 
FEMALES MALES NUMBER 

PARENTS BENT-TAIL NORMAI BHENT-TAII NORMAL OF LITTERS 
bn bn 2 *& 

Bent-tail o% 69 56 14 
Bnbn 2 

Bent-tail 29 
Bn Bn G 


? 
normal 
Bn Bn 2 &X 
Bent-tail 
Bn bn 2 X 
normal 29 45 37 18 


A 


Heterozygous females did not appear to differ from normal females as 
regards rate of growth; approximately 50°) of the Bent-tail males and all 
the homozygous females were smaller at birth and matured more slowly 
than normal mice. Four Bent-tail males died within a few months; the 
toes of the forefeet of these males had been turned against the foot causing 
the mice to hobble on their “knuckles.”’ 

Genetic Analysis.—-Matings between normal females and Bent-tail 
males yielded litters containing only Bent-tail females and normal males 
(mating |, table 2); matings between heterozygous females and Bent-tail 
males produced litters containing normal and Bent-tail males and only 
Bent-tail females (mating 2, table 2). These observations indicate that 
Bent-tail is a dominant, sex-linked mutation. The results of these matings 
also suggested that the parental Bent-tail males were either heterozygous 
(Bn bn) or hemizygous (Bn 0) with Bn on the X-chromosome. 

Matings between heterozygous females and Bent-tail males were ex 
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pected to yield 2 types of Bent-tail females: homozygous (Bn Bn) and 
heterozygous (Bn bn). Four females with tails half normal length and 
similar in appearance to the tails of Bent-tail males were found among the 
first 19 Bent-tail females from matings between heterozygous females and 
Bent-tail males and were found to be homozygous on the basis of the 
following observations. Two of the 4 females were fertile, yielding only 
Bent-tail progeny when mated with either normal or Bent-tail males 
(matings 3 and 4, table 2); matings between these females and Hent-tail 
males produced only Bent-tail females with tails half normal length. 

Bent-tail males mated with C3H females, kindly furnished by Dr. K. B. 
DeOme, University of California, yielded only Bent-tail females and nor- 
mal males indicating that the expression and inheritance of Bent-tail were 
not influenced by the C3H genotype. 

A significant deviation between the observed and theoretical ratios was 
noted in litters from matings between heterozygous females and Bent-tail 
males (mating 2, table 2); the normal class comprised approximately 50°F 
instead of 25°, of the population. Since Bent-tail males and homozygous 
Bent-tail females do not appear to be as vigorous as either normal males 
and females or heterozygous females, the deviation between the observed 
and theoretical ratios may be due to a lethal effect of Bent-tail in the male 
or homozygous female. 

The mean litter size was determined for different types of matings in- 
volving Bent-tail (table 2): (1) bn bn 2 X Bent-tail o—8.9; (2) Bn bn 
9 & Bent-tail o& -6.3; (3) Bn Bn? XK normalc’ -4.3; (4) Bn Bn? X Bent- 
tail o'-4.5; (5) Bn bn 2 X normal o -7.8. These data suggest that hetero- 
zygous females (Bn bn) are less fertile than normal females and homozy- 
gous females (Bn Bn) are less fertile than heterozygous females. 

Discussion. Sex-linked mutations may be distinguished from autosomal 
mutations by their association with sex in certain types of matings. If 
the mutation affects a visible character, the typical crisscross inheritance 
is Observed. If the mutation manifests itself only as a lethal for males, a 


discrepancy in the sex ratio is noted.'! A third method has been considered 
by Fisher, ef a/.6 who interpret recombination values for sex and certain 
mutations in mice as evidence for partial sex-linkage. Bent-tail is in 
herited in the classical pattern of sex-linkage found in Drosophila melano 


gaster. 

Does the Y-chromosome in the house mouse carry an allel for Bent-tail? 
Since Bent-tail is dominant, it would be difficult to distinguish between a 
recessive allel (heterozygosity) or no allel (hemizygosity) in the Y-chromo- 
some. If Bent-tail were in the homologous segment of the X-chromosome 
and a recessive allel in the Y-chromosome, crossing over in such a male 
would yield a Y-chromosome with the Bent-tail gene and would be de- 
tected in the progeny of matings between normal females and Bent-tail 





Vo. 38, 1952 GENETICS: E. D. GARBER 879 


males by the presence of a Bent-tail male. Matings designed to test this 
possibility have not yet yielded a Bent-tail male. These negative results 
may be due either to an extremely low frequency of crossing over or to 
hemizygosity with the Bent-tail locus in the non-homologous segment of 
the X-chromosome. 

The expression of Bent-tail in males and females furnishes evidence sup- 
porting the hypothesis that Bent-tail males are hemizygous for the Bn 
locus. Although heterozygous females have tails normal in length, both 
homozygous females and Bent-tail males have tails half the normal length. 
If the expression of Bent-tail were dependent on a dosage effect, heterozy- 
gosity and homozygosity would be expected to yield dissimilar phenotypes. 
Since Bent-tail males are similar in phenotype to homozygous Bent-tail 
females and much more extreme than heterozygous Bent-tail females, 


hemizygosity of the Bent-tail males is indicated. Until such time as 
genetic data may be obtained indicating that crossing over has introduced 
the Bn locus into the Y-chromosome, Bent-tail males will be considered to 
be hemizygous (Bn Q) for this mutation which is probably located on the 


non-homologous segment of the X-chromosome. 

Summary.—-Bent-tail (Bn) is a dominant, sex-linked mutation in the 
house mouse visibly affecting the structure and length of the tail as well as 
unknown characters leading to inviability in males and homozygous fe- 
males. Penetrance is high, at least 95°; in heterozygous females and 100% 
in males and homozygous females. Homozygous temales and Bent-tail 
males have tails half normalin length. Bent-tail males are considered to be 
hemizygous (Bn 0) for the Bn locus which may be in the non-homologous 
segment of the X-chromosome. 
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THE INSOLUBILITY OF SETS OF DIOPHANTINE EQUATIONS 
IN THE RATIONAL NUMBERS 


By Nesmitu C, ANKENY 
MATHEMATICS DEPARTMENT, JOHNS HOPKINS UNIVERSITY 
Communicated by S. Bochner, August 8, 1952 


I 


Denote by M(7, a, a2, ..., @,) = M(7) the set of all positive integers 
m < 7 for which the equation 


a,X ,™ + A,X 2™ ‘ +. riage = () (1) 


has a non-trivial solution in the rational integers. By a non-trivial solu- 
tion we mean a solution of (1) where not all of the X's are zero. 

Putting an assumption on qd, ..., @, we can deduce an upper bound on 
the density of A/(7). 

Assume 


n 
2 8, * 0, (2) 
1 

where &, +1 or o. 

TueoreM 1. If a, ... , dy satisfies (2), then M(T) = o(T), as T > &, 
t.e., for almost all m equation (1) has no non-trivial solution. 

One rather important case that (2) excludes is when n = 8, a) = dy = 
d; = 1. However, some information can be gathered but only in the “‘first 
case’ of Fermat’s theorem. 

THroreMm 1], The density of integers m < 7 for which the equation 


Xi" + X," + X3" = 0 


has a rational solution and for which (X,XoX3,m) = 11s 0(T) as T > @. 
One can extend these results from the coefficients a), ..., @, contained in 
the ground field to similar results about a, , a, contained in any alge- 


braic number field /. However, in this case one must change condition (2) 
to 


n 


>» ar; ¥ O, 
j= | 
where A, = O or any root of unity contained in F. 

In this present paper the complete proof will not be given. We shall 
give the main supporting Lemmas, a sketch of their proof, and in conclu- 
sion, the deductions of Theorems I and II from the Lemmas. The Lem- 
mas are of two different natures, namely from algebraic number theory and 
from analytic number theory. 
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Notations: (,, C2, ..., €z will be positive constants. 7), U will be large 
positive variables tending to + ©. p will denote all primes #1 (mod n!). 
g, m, r and ¢ denote positive integers. 


I] 


LemMA 1. If a, ..., Q, satisfies (2), g > 2, and (g, n!) is etther 1 or . 
then ¢, cannot satisfy an equation whose coefficients are ay, dz, ..., An. 

Proof: The proof consists in induction on the number of prime factors 
of g. If gisa prime the Lemma is immediate. Proceed by induction, but 
transfer the ground field R to other cyclotomic fields. 

LEMMA 2. If a, dz, ..., ad, satisfy condition (2), then for a given m 
there exists no non-trivial rational solutions of (1) provided we can find a 
rational prime p such that 


m divides p — 1, mr = p — 1, (3) 
(y, ni} = 1 or 2, (4) 
g(r) << a 1 log p, (5) 


where a = log (Ja,| + |d2) +... + |a,|), and g(r) is the Euler ¢ fune- 
tion. 

Coroiiary. If n = 3,4; = dy = a3 = 1, then (1) has no non-trivial solu- 
tion relatively prime to m (see Theorem II for the ‘first case’) provided a prime 
p exists that satisfies (3), (4) and (3). 

M. Krasner has independently proved the Corollary. His proof is most 
elegant and elementary. A less general form of Lemma 2 was proved by 
H. S. Vandiver. 

The general proof proceeds as follows: Assume there exists a p which 
satisfies (3), (4) and (5), and that (1) has a rational solution such that 
X,X_... X, 4 0 (mod p). Without loss of generality, assume (X,, Xo, 

., X,) = 1. Then consider (1) in the field R(¢,). 

As p = 1 (mod r), the ideal factorization of p is (p) = P,P, ... Py, in 
R(¢,), where s = g(r), and Naw), (P;) = p. Hence, the group of mth 


power residues of the multiplicative cyclic group of residues (mod P,) has 
(p — 1)/m = relements. One shows that the elements ¢//, 7 = 0,1, ..., 
r — | are incongruent (mod P,). So ¢/ form a subgroup of r elements in 


multiplicative subgroup of residues (mod. P,). Hence, these two sub- 
groups must coincide. 


As 
aX,” + + a = # m= Q, 
a fortiori, 


aX," Si ree + re == (0 (mod FP); 
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or by the coinciding of the two subgroups, 
ag," + + a,¢," Q (mod P;) 


for some n-tuple of integers (4, , tn). Hence, 


p = Nui) n(P)) 


i as , 5 @ > oe 
divides Nr ¢,).a(aig,' + Gt; ). ret, 


r ge te bn 
Nace) nlat,' + + anf,") S ({ai| + + 


pS (\a,) + + |a,|)*”" 


which is a contradiction to our hypothesis, unless 
w t « tn 
aig, + + nhs Q). 


This case, however, is impossible by Lemma 1. 
So we have shown that 


X1Xe ea 0 (mod p). 


Hence, p divides one of the variables, say X,. However, proceeding in the 
same way with the truncated equation a,X," + a2Xo" + ... + dyiXy,-1" 
we will see that p will divide each X;, 7 = 1, 2, ,n. This is a contra- 
diction to (X;, X», ,X,) = 1. This proves Lemma 2. 

In the case when n = 3, ad; = dy = dy = 1, then the above proof shows 
that p will divide X,y. If (m, XiX2X;) 1, we have by using the well- 


known criterion on Fermat’s equation, that 

pt! (mod gq’) 
for all primes g dividing m. As p | (mod m), this implies p == 1 (mod 
m*). But 

p —1=mr<ma' log p S m log p. 
So m > p/log p, or p < 2m log m, or p < m?*, p # 1 (mod m?) a contradic 
tion, proving the Corollary. 
II] 


The following Lemmas use the tools of analytic number theory, namely 
the Brun-Selberg sieve methods and certain “probability” theorems re 
garding the density of numbers having a fixed number of prime factors. 

LemMMA 3. The density of primes p satisfying: 


cr Ul <p < 2arl, (6) 
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where L’ is fixed and r varies over all positive integers < log U, 
p | (mod r), (7) 


all prime factors of p — 1/r are greater than UL’, ts greater than c,U'/log U. 

Proof: The method of proof of Lemma 3 is essentially the same as used 
in Goldbach’s conjecture, showing that every even number is the sum of a 
prime plus a number not containing more than |1/¢;| prime factors. 

In Estermann’s' paper the constant c; can be replaced by c; 6, if 
one assumes the Extended Riemann Hypothesis. The recent advances 
of Linnik on both Goldbach’s conjecture and primes in arithmetic progres 
sion make the above assumption superfluous, but one must replace 1,6 by 
an unknown constant ¢;. 

LemMaA 4: The density of primes p satisfying: 


p < 2erU, (S) 
| (mod g) (9) 


where g varies over all numbers all of whose prime factors are > LU and g lies 
in the interval o,U <g <cgl’; rvaries from 1 to log U, is less than |c7(cy — 65) 
U|/log U. 

Proof: The proof of Lemma 4 is relatively easier than Lemma 3 as it 
only demands an upper bound on a set of integers on which various multi 
plicative conditions have been imposed. 

From Lemmas 3 and 4, one deduces that a fixed finite proportion of all 
numbers g which satisfy condition (4) of Lemma 4 also satisfies the condi 
tions imposed on m in Lemma 2 with 2c. a '. Hence, for this finite 
proportion of g we see that equation (1) has no non-trivial solution. We 
add to this the obvious fact that 1f (1) has no solution for m, (1) has no 
solution for any multiple of m. 

So we now have reduced the proof of Theorem I to showing that almost 
all numbers m < 7 are divisible by a fixed finite proportion of all numbers 
g which satisfy (4) of Lemma 4 where l’ varies from 100 to 7. 

If c; 2 1/2, we would then see that g would necessarily be a prime. 
Hence, our problem would be reduced to finding the density of numbers 
which are not divisible by a set of finite density of primes. This can easily 
be done by the sieve method. If ¢; is just a given positive constant the 


problem can be solved by combining the sieve method by certain “‘probabil- 


ity arguments” in number theory. We show that if a number ¢ has more 
than (log log ¢)’’'* prime factors, then ¢ has probability | of being divisible 
by one of our Set g. Then by use of “probability arguments’? we can show 
that the numbers which have fewer than (log log ¢)"’ prime factors have 
density zero. 

This completes the proof of Theorems I and II. 
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When R is replaced by any algebraic number field F then almost all argu- 
ments of the above proof carry over. The rational primes p are replaced 
by the primes of first degree in F (i.e., the primes which completely split in 


F). 


! Estermann, T., ‘Eine neue Darstellung und neue Anwendungen der Viggo Brunchen 
Methode,” J. Reine Angew. Math., 168, 107-116 (1932), 


ON THE NORMAL ORDER OF ARITHMETIC FUNCTIONS 
By RICHARD BELLMAN AND HAROLD N. SHAPIRO 


THE RAND CorporaAtION AND DEPARTMENT OF MATHEMATICS, 
NEW YorK UNIVERSITY 


Communicated by H. S. Vandiver, August 20, 1952 


1. Introduction, -Given an arithmetic function /(”) possessing a mean 
value in the sense that }> /(m) ~ ex, ¢ ¥ 0, asx — ©, the question was 


nox 
raised by Gauss as to the persistence of the mean value if one takes the 
moving average > f(n)/y, where y = y(x) + © asx — ©. Two 
xSnsery 
interesting problems of this genre are those of the distribution of squarefree 
numbers in small intervals, and, in particular, the determination of the 
inaximum distance between consecutive squarefree numbers, and the similar 
question for numbers which are the sum of two squares. The first problem 


. . . - ie » ® 4 
involves the estimation of the sum os u*(n) while the second re- 


£2RKAR2 + F 
quires the estimation of > rin). 
s-aSsaery 
More generally, we may take a function such as d(m), the Dirichlet di- 
visor function, which has the mean-value log v and seek to determine the 
functions y(x) for which > d(n)~ ylogx asx —> &. It is clear 
xansxty 


that any error term in the mean-value formula such as }> f(m) = cx + 
“aa 


O(x*), 0 << a@ < 1, permits one to assert that 3. [(n) ~ cx’ for b> a. 


tee ee 
In many cases, as in the distribution of primes and squarefrees, one does 
not want the full force of the asymptotic relation, but merely a relation of 
the type :> f(n) > 0. Thus, for example, Hoheisel proved there 
ee. Ie ° 
was a prime between x and x + x%,a = 1-! 33,000, for large x, a result which 
cannot be derived by use of the present estimates for the remainder in the 
prime number theorem. This result was improved considerably by In- 
gham!' and his results were applied to investigation of moving averages of 
arithmetic function by Fogels. Unfortunately, in general, the compli- 
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cated analytic techniques do not yield results anywhere near as precise as 


those derived by more elementary means. 
The purpose of this note is to indicate a technique for deriving the normal 
order of F(x, y) = > f{(n) for a large class of arithmetical func- 
SSNS T 9 
tions. The normal order of F(x, y) is defined to be g(x) if for any « > 0, 


“ 


we have (1 — e)g(x) < F(x, vy) < (1 + e)g(x) for all x < 2 except for a set 


of total number o(z) asz—~> ©. 

For the case of squarefree numbers we obtain a particularly striking re 
sult which is stated below. 

2. The Method..-The technique we employ is intimately connected to 
the methods introduced by Schnirelman*? and Turan.’ We consider the 
sums 


x * 


x < 


ys f(n)), 
r+y 


n-< 


bi fin)),’ (1) 
r+y 


% 5 


where y(x) is a monotone function of x which +> © asx-—» @. It is easy 
to verify that S\(z) ~ }> y(x) if } f(n) ~ x. The evaluation of S:(z) is, 
x>@ ns 
however, more difficult, involving the asymptotic behavior of the sums 
T(z, k) = > f(n) f(n + k). This type of sum was first encountered by 
n—& 

Ingham‘ in connection with the question of the mean-value of ¢(s) on the 
critical line, with f(m) = d(n). 

If we prove that S.(z) ~ 3S y?(x), it will follow, with some slight addi- 


tional assumptions concerning y(x), that F(x, y) has the normal order y. 
In this way we can show that 7 r(n) has the normal order 4y/ 2? 


S>a>RaAA2T ¥ 
if y = (log x)’, and that = d(n) has normal order log x if y = 
£23822 7T 9 
(log x)'**, a > 0. Interestingly enough, the result is not true if y = log x. 

The proofs of these and many further results will be presented in a sub- 
sequent publication. 

3. The Squarefree Problem. For the case of the distribution of square 
free integers we can obtain a very sharp result using only the evaluation of 
Si(z) and the following lemma, itself of interest. 

LEMMA. Guivenany «> 0, we have 


u*(n) < (1 4+ 
for x > xo(e), where y is any function of x which approaches ~ asx— @, 
This is a simple consequence of the sieve method of Brun. Combining 
this with the result 
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le» 


O25 
p(n))~ — d& y(x), 
y Wyse 


we have 


THEOREM |. /f y(x) 1s any function such that lim y(x) = ©&, then 


eee 
the number of n < z such that the interval (n, n + y(n)) contains no square- 
free integer is O(2) AS Z-—> &, 


THEOREM 2. /f lim y(x) wo, and for some X,0 <» <1, y(x)/y(Ax) 


—- 


is bounded, then the normal order of ‘ u(n) is (6/m*)y. 


x noxtry 


It is interesting to compare these results with the results which hold for 


‘ Ox ele 
all nm. ‘The formula : > p(n) , + O(Y x) furnishes [x, x + kv/ x} as 
BS x Tv 


an interval which always possesses a squarefree number for & sufficiently 
large. It was proved independently by a number of mathematicians that 
this could be whittled down to [x, « + kx'’'|, and even to [x, x + kx'/*/ 
log x}. Recently, by means of a novel and ingenious method, Roth® 1m- 
proved this to [x, x + x" * ‘|. Although the result is undoubtedly true 
for [x, x + x], the result seems very difficult to establish. In the other 
direction it is easy to show that there are infinitely many intervals [x, x + 
¢ log x/log log x| which do not possess a squarefree number. 

' Ingham, A. E., Quart. J. Math., 8, 255-266 (1987). 

2 See, for example, E. Landau, l’ber einige neuere Fortschritte der additiven Zahlen- 
theorie, Cambridge, 1937 

§ Turan, P., J. London Math. Soc., 9, 274-276 (1934). 


‘Ingham, A. E., J. London Math, Soc., 2, 202-205 (1927) 
» Roth, K. F., Lbid., 26, 263-268 (1951) 


A SOCIAL EQUILIBRIUM EXISTENCE THEOREM* 
By GERARD DEBREU 
CowLks COMMISSION FOR RESEARCH IN ECONOMICS 
Communicated by J, von Neumann, August 1, 1952 


In a wide class of social systems each agent has a range of actions among 
which he selects one. His choice is not, however, entirely free and the ac- 
tions of all the other agents determine the subset to which his selection is 
restricted. Once the action of every agent is given, the outcome of the so- 
cial activity is known. The preferences of each agent yield his complete 
ordering of the outcomes and each one of them tries by choosing his action 
in his restricting subset to bring about the best outcome according to his 
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own preferences. The existence theorem presented here gives general con 
ditions under which there is for such a social system an equilibrium, 1.e., 
a situation where the action of every agent belongs to his restricting subset 
and no agent has incentive to choose another action. 

This theorem has been used by Arrow and Debreu* to prove the existence 
of an equilibrium for a classical competitive economic system, it contains 
the existence of an equilibrium point for an N-person game (see Nash* and 
Section 4) and, naturally, as a still more particular case the existence of a 
solution for a zero-sum two-person game (see von Neumann and Morgen 
stern, Ref. 11, Section 17.6). 

In Section | the topological concepts to be used are defined. In Section 
2 an abstract definition of equilibrium is presented with a proof of the 
theorem. In Section 3 saddle points are presented as particular cases of 
equilibrium points and in connection with the closely related MinMax 
operator. Section 4 concludes with a short historical survey of results 
about saddle points, fixed points for multi-valued transformations and 
equilibrium points. 

Only subsets of finite Euclidean spaces will be considered here. 

1. Topological Concebts..Two sets in R” are said to be homeomor phi 
when it is possible to set up between them a one-to-one bicontinuous (/ and 
h~' continuous) correspondence / (called a homeomorphism). 

A convex cell C in R" is determined by r points 2*(k l, ,r); itis the 
set 

C = {siz pd ye oc, = Ofork a 2 >> re 
k=! 1 
Such a set is closed. 

The product of two convex cells A ¢ R' and B ¢ R” is @ convex cell 
Cc R't™. Let A be generated by the p points x! (7 ¥. , Pp) and B 
by the g points y’ (7 = 1, ,g). Denote by C the convex cell in R'*” gen 
erated by the pq points (x‘, vy’). Obviously A &K B 2 C, and one shows 
easily that A KX Bc C. 

A geometric polyhedron is the union of a finite number of convex cells in 
R". It is clearly closed. 

The product of two geometric polyhedra P, Q is a geometric polyhedron. 

4 qd 


Let P = / yu B ) where the A, (the B,) are convex cells in 
1 


) 


R' (in RR"). The relation P XK Q (uw AD) xX Cae) (Art B) 


uw 


proves the result. 

A polyhedron is a set in R" homeomorphie to a geometric polyhedron 
(called geometric antecedent of the first one). 

The product of two polyhedra is a polyhedron (since 1t is homeomorphic 
to the product of the two geometric antecedents). 
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Let J = {t}0 < t S 1} denote the closed interval [0, 1] on the real line. 
A nonempty set Z of R" is said to be contractible, or more precisely, deform- 
able into a point 2° ¢ Z, if there exists a continuous function //(t, z) (called 
a deformation) taking J X Z into Z such that for all z « Z, 7(0, z) = 2 and 
H(i, zs) = 2°. 

The product of two sets X ¢ R‘, Y ¢ R”™ deformable into the two points 
x° « X, y° « Y, respectively, is clearly deformable into the point (x°, y°). 


0 


Finally the real function t = ~~ 
é 


increasing from —1 to +1 when 0.. -reases from — © to +o. It estab- 
lishes a one-to-one correspondence bev.veen the closed interval [—1, +1] 
and the set R of all real numbers to which are added two elements — © and 
+. Open sets in R are defined as images of the usual open sets in [—1, 
+1], an order is defined in R as an image of the usual order in [—1, +1]. 
R endowed with this topology and this order is called the completed real line 
(which can naturally be defined directly®). 

2. Equilibrium Points.-Let there be v agents characterized by a sub- 


of the real variable 0 is monotonically 
: ) 


0 a}. 


scripte = 1, ..., v. 

The «th agent chooses an action a, in a set Y,. The v-tuple of actions 
(a,, ..., a@,), denoted by a, is an element of % = %X... XY. The payoff 
to the cth agent is a function f(a) from Y% to the completed real line. 

Denote further by d, the (v — 1)-tuple (a), ...,@,~1, @,41, ...,a,) and by 
Wf, the product WX... KM) K Wa &... KY, Given a, (the actions of 
all the others), the choice of the cth agent is restricted to a non-empty, com- 
pact set A(a,) ¢ YX; the eth agent chooses a, in A,(d@,) so as to maximize 
f(a, a,), assumed to be continuous in a, on A,(d,). 

This background makes the following formal definition intuitive : 

Definition a* is an equilibrium point if for all « = 1, ..., va. € A,(a°) 
and f(a*) = Max f,(a;, a). 

au € Au (ae*) 

The graph of the function A((d,) is defined as the subset of %, X Y, 
G, = |(@, a,)\a,¢€A(a)}. For any a, A,(a,) is always understood to be 
non-void. 

Turorem. Foralli = 1, , v, let UX, be a contractible polyhedron, A (a,) 
a multi-valued function from W, to 1, whose graph G, ts closed, f,a continuous 
function from G, to the completed real line such that ¢f{a,) = Max 


ae Ag (ai) 
f(a, a,) ts continuous. If for every « and a, the set Ma, = {a, € A,(a,) 
f(a, a,) ¢,(4,)|} ts contractible, then there exists an equilibrium point. 

The proof uses as a lemma a particular case of the fixed point theorem 
of S. Eilenberg and D. Montgomery® or of the even more general result of 
E. G. Begle.* 

Let Z be a set and ¢ a function associating with each z « Z a subset $(z) 
of Z. We have defined above the graph of ¢ as the subset of ZX Z, | (z, z’)| 
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2’ep(z)}. is said to be semicontinuous if its graph is closed. <A fixed 
point of dis a point 2* such that 2* € ¢(2*). 

LemMa. Let Z bea contractible polyhedron and @: Z—» Za semicontinu- 
ous multi-valued function such that for every z € Z the set o(2) is contractible. 
Then o has a fixed point.t 

Yl, the product of v contractible polyhedra, is a contractible polyhedron 
(Section 1). Define on Y% the multi-valued function @¢ as follows: 


¢(a) = M, xX... X Ma, 


Since \/,, is contractible for all « and a,, ¢(a) is contractible for all a € Y 
(Section 1). To be able to apply the lemma it remains only to show that 
¢@ is semicontinuous. 

For this first define in {, X Y%, the set 


M, = {(d,a,)|a, € M,,}. 


‘ 


The equivalent definition 


M - ;(a,, a 


‘ 


J €G f(a, a) ¢,(d,)} 


shows that 17, is closed since G, is closed and f, and ¢, are continuous. 
The graph I of ¢ is the subset of YU X YF 


r = {(a,a’)\a’ € d(a)} = {(a,a’) a,e M,, for alle} 
f(a, a’)| (a, a.) € M, for all c}. 
Consider the subset of % & YF 
M, = {(a, a’)| (4, a,) « M,}; 


WW, is closed since M/7,is. AsT = a Wi, T is closed. 


The conclusion of the lemma is then that there exists a* ¢ Y{ such that 
a* ¢ p(a*), ie., for alle, a” « Mat; this is the definition of an equilibrium 


point a*. 


The requirement that ¢,(@,) be continuous is a joint requirement on the 
two functions f, and A,(@,);_ it is therefore not well adapted to applica- 
tions. The following Remark tries to overcome this. 

The function A,(d,) is said to be continuous at a if for any a? e A,(a’) 
and any sequence (d@") converging to a’, there exists a sequence (a/) con 
verging to a such that for all n, a" € A,(@)). 

Remark: If A<a@,) has a compact graph G, and is continuous at a), if f, 
is a continuous function from G, to the completed real line, then ¢(a,) is con- 
tinuous at a. 

We drop subscripts « everywhere and reason as if f took its values in the 


; eid I 
real line (the isomorphism between the completed real line and the 
f 
er] 
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closed interval |—1, +1] immediately extends the results to the general 
case ). 

(a)  L’sing only the compactness of G and the continuity of f we first prove: 
For any sequence (d") converging to a’ and any e > O, there is an N such 
that n > Nimplies g(a") < g(a°) + e€ (in other words, g(a) is upper semi- 
continuous at a”), 

For every n, choose a" « A(a") such that f(a", a") = g(a"). Since G 
is compact it is possible to extract from the sequence (a”, a") a subsequence 
(a”"’, a"’) converging to (a°, a°). 

By the continuity of /, f(a’’, a’) |which is = g(a"’)| tends to f(a", a°) 
|which is S g(a@°)|. Therefore there exists N’ such that n’ > N’ implies 
¢(a"’) < g(a") + «. Since from any sequence (a”") converging to @° it 1s 


possible to extract a subsequence (a@” ) having the desired property, any se 


quence (a") converging to a° has the property. 

(8) UL’sing in addition the continuity of A(a) at a° we prove: For any se 
quence (d") converging fo d° and any e > 0, there is an NV such that n > N 
implies g(a") > ¢(a°) « (in other words, g(@) is lower semicontinuous 
at a”). 

Choose a° ¢€ A(a°) such that f(a", a) ¢g(a°). By continuity of A(a) 
at a", there is a sequence (a") converging to a° such that for all n, a” 
A(a"). By the continuity of f, f(a", a") [which is < g(a")! tends to 
f(a", a°) {which 1s = g(a°)|. Therefore there exists NV such that n > N 
implies g(a") > g(a°) — «€. 

(a) and (8) together naturally prove that g(a) is continuous at @°. 

3. Saddle Points and MinMax Operator. In this section x « X ¢ R‘, 
ye Y ¢ Rand f(x, y) isa function from X X& Y to the completed real line. 
A saddle point of fis a point (x°, y°) such that 


Min f(x°, y) ier, 9") Max f(x, y°). (1) 


y x 


It is a very particular case of an equilibrium point for two agents: 


em YI A A (dy) x. f(a) == f(x, y) 
ay y, YW, =: ¥, Ao(d») r, fo(a) = —f(x, y) 


One obtains therefore by using the remark: 
Coro._iary. Let X, ¥ be two contractible polyhedra, and f(x, y) a con 
tinuous function from X X Y to the completed real line. If for every x°« X, 


Uys fy o \ f(e%, vy) = Min f(x°, y)} ts contractible and for every y® « Y, 
ye V 


Vyo = fxeX) f(x, y°) = Max f(x, v°){ 7s contractible, then f has a saddle point. 
re X 


This corollary contains as more and more particular cases the saddle 
point theorems of Kakutani,’ von Neumann (Ref. 9, p. 807, and 10), and 
von Neumann and Morgenstern (Ref. 11, Section 17.6). 
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The special interest of saddle points comes from their intimate relation 
with the Min Max operator. 


From now on X, Y are assumed to be compact and f{(x, y) to be continuous. 
We know from the Remark that Min /(x, y) [resp., Max /(x, y)] 1s@ con 
y x 

tinuous function of x {resp., vy}. The following results, already given in 
Ref. 11, Section 13, are proved here for completeness. 

(a) Max Min f(x, y) S MinMax f(x, y). 

x y y x 
Let A = {x’| Min f(x’, y) Max Min f(x, y){, B fy’! Max f(x, y’) = 


x x 


Min Max /(x, y){. 


y x 

Iix’ «A and y’ € B, 
MaxMin f(x, y) = Min f(x’, y) S f(x’, y’) S Max f(x, y’) 
MinMax/(x, y). 


\ 


The result follows from a comparison of the first and last terms. 
(b) The existence of a saddle point (x°, y°) implies the equality 


Max Min / Min Max /| 


x y x 


From the definition | it follows that 
Max Min f(x, y) 2 Min f(x°, y) (x, 3°) 


Max f(x, y’) 2 MinMax f(x, y) (3) 


x V x 


which together with (a) gives the result. It also gives Max /(x, y") 
MinMax /(x, y) i.e., y® « B, and similarly x° € A. 
y 3 
(ce) The equality MaxMin f = MinMax f implies the existence of a saddle 
point. al na 
Assume that the equality holds and take x° « A, y® e B, Eq. (2) gives Min 
. ¥y 


y) = f(x®, y°) = Max f(x, y"), which is the definition Eq. (1) of a saddle 


(x? 


, 
point (x°, y°). We have, incidentally, proved 

(d) the set of saddle points is either empty or equal to A X B. 

4, ILlistorical Note. A function f(z) from a set Z to the completed real 
line R is said to be quasi-convex (resp., guasi-concave) if for any a ¢ R, the 


set of z eZ such that f(z) S a (resp., f(z) 2 a) is convex. 


" 


Let S, = is « R"\s, 2 Ofork = 1, , nm and 
k 
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In his first study on the theory of games, J. von Neumann’ proved: 

(1). Let f(x, y) be a continuous real-valued function for x ¢ S, and y € S,,. 
If for every x° € S, the function f(x°, y) is quasi-convex, and 1f for every y° 
eS, the function f(x, y°) 1s quasi-concave, then f has a saddle point. 

In another paper on economics" he later proved a closely related lemma 
which S. Kakutani’ restated in the more convenient form of the following 
(equivalent) fixed point theorem: 

(II). Let Z bea compact convex set in R" and @: Z—> Za semicontinuous 
multi-valued function such that for every z € Z the set $(z) 1s non-empty and 
convex. Then has a fixed point. 

The convexity assumptions were, however, irrelevant and S. Eilenberg 
and D. Montgomery® gave a fixed point theorem where convexity was re- 
placed by acyclicity. Their result was further generalized by E. G. 
Begle.’ 

These last two theorems deserve particular attention as valuable con- 
tributions to topology whose origin can be traced directly to economics. 


The notion of an equilibrium point was first formalized by J. F. Nash* 
in the following game context. There are v players; the ith player chooses 
a strategy s,in S,,; his payoff is f,(s), a polylinear function of s;, ..., 5,. 
An equilibrium point is a y-tuple s* such that for all c, f,(s") = Max 


s S 
‘ e nm, 


ci, 4 Nash proved the existence of such an equilibrium point. 


iJ* 


* Based on two Cowles Commission Discussion Papers, Mathematics 412 (Nov. 1, 
1951) and Economics 2032 (Feb. 11, 1952). This paper has been undertaken as part 
of the project on the theory of allocation of resources conducted by the Cowles Com- 
mission for Research in Economies under contract with The RAND Corporation. To 
be reprinted as Cowles Commission Paper, New Series, No. 64 

One of the main motivations for this article has been to lay the mathematical founda 
tions for the paper by Arrow and Debreu;? in this respect I am greatly indebted to K. J. 
Arrow. Acknowledgment is also due to staff members and guests of the Cowles Com- 
mission and very particularly to [. N. Herstein and J. Milnor. I owe to J. L. Koszul 
and DD. Montgomery references 6 and 38. Finally I had the privilege of consulting 
with S. MacLane and A. Weil on the contents of Ref. 6. 

t The statement of E. G. Begle (Ref. 3, p. 546) is indeed much more general and the 
existence theorem can accordingly be generalized. Instead of a contractible polyhedron 
one might take for example an Absolute Retract (as defined in |Ref. 4, p, 222]) using the 
fact that the product of two A.R. is an A.R. !Ref. 1, p. 197]. For finite dimensions 
“Absolute Retract”’ is equivalent to ‘contractible and locally contractible (Ref. 4, pp. 
235-236) compact metric space’ [Ref. 4, p. 240] 

' Aronszajn, N., and Borsuk, kK., “Sur la somme et le produit combinatoire des ré 
tractes absolus,’’ Fundamenta Mathematicae, 18, 193-197 (1982) 

* Arrow, K. J., and Debreu, G., ‘Existence of an Equilibrium for a Competitive 
Economy,'’ Econometrica, in press (1953) 

§’ Begle, E. G., “A Fixed Point Theorem,” Ann, Math., 51, No. 3, 544-550 (May, 
1950) 

4 Borsuk, K., ‘Uber eine Klasse von lokal zusammenhiangenden Raumen,” Funda 
menta Mathematicae, Vol. 19 (1932), p. 220-242. 





VoL. 38, 1952 MATHEMATICS: E. G. KUNDERT 893 


5 Bourbaki, N., Eléments de Mathématique, Premiére partie, Livre III, Chap. IV, 
§4, Hermann, Paris, 1942 

6 EKilenberg, S., and Montgomery, D., ‘Fixed Point Theorems for Multi-valued 
Transformations,’’ Am. J. Math., 68, 214-222 (1946 

7 Kakutani, S., ‘‘A Generalization of Brouwer’s Fixed Point Theorem,’’ Duke Math. J., 
8, No. 8, 457-459 (September, 1941) 

8 Nash, John F., “Equilibrium Points in N-Person Games,’’ Proc. Nati. ACAD. 


Scr. 36, 48-49 (1950) 

® Neumann, J. von, ‘Zur Theorie der Gesellschaftsspiele,”’ Math. Ann., 100, 205-320 
(1928) 

10 Neumann, J. von, “Uber ein 6konomisches Gleichungssystem und eine Verall 
gemeinerung des Brouwerschen Fixpunktsatzes,”” Ergebnisse eines Mathematischen 
Kolloquiums, 8, 73-83 (1937), (translated in Rev. Economic Studies, XII, No. 33, 1e9 
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A RELATION BETWEEN POLES AND ZEROS OF A SIMPLE 
MEROMORPHIC DIFFERENTIAL FORM AND A CALCULATION 
OF CHrRN'S CHARACTERISTIC CLASSES OF AN ALGEBRAIC 
VARIETY 
By E. G. KUNDERT 
DEPARTMENT of MATHEMATICS, LOUISIANA STATE UNIVERSITY, BATON ROUGE, LA. 
Communicated by S. Lefschetz, August 18, 1952 

It is well known! that on a Riemann surface the degree of the divisor of a 
differential is always = 2p — 2 = x. We give here a generalization of 
this formula for simple meromorphic differential forms of higher dimensions. 

Let @ be such a differential form on the complex manifold .W". We de 
note by c(a) its zero cycle and by s(a@/ its polar cycle, and by C(a@) and 
S(a) the cohomology classes of the duals of cla) and s(a). 

THEOREM I: C(a@) can be expressed as a polynomial of degree n in S(a) in 
the cohomology ring of AM": 


k 


n 
Cl ap = Aad a) + z; rl, k+1° Sa) 
k=] 


where the cocycles , are the covariant Chern classes of M". 

Not only because of Theorem I, but for various other purposes, it 1s de- 
sirable to have a simple method to calculate the Chern classes. We give 
such a method for the case of an algebraic variety. 

THEOREM II: The Chern classes of an algebraic variety can be determined 
by the following formulae: 
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n 
( 1)*-T, = i," 4 2. IT .—241° TT,” x 
k=l 


o%- YT (I,) 
1)! Oll,, 


where \\, denotes the cohomology class determined by the dual to the intersection 


cycle P* \ V" in the projective imbedding space P”*'; P* being a k-dim. 
(k 1) 
plane in P"*', The symbol in (3) has to be interpreted as the formal 


0” 
ke — 1" derivative of the polynomial (2) with respect to Ip. 

Remark: One observes that (2) implies a calculation of the Euler char 
acteristic of 1". 

Proof of Theorem 1: Tf ais the form ya dz, and a, is locally written 
as the quotient p;/q, then denote by qg the lowest common multiple of the 
qi. We attach to the form a the field of covariant line-elements F defined 
by }q, Pu » Pn}. Cla) is then-—as one at once confirms—nothing 
else than the obstruction of the field F. We introduce now as a reference 
field, the field of line-elements G corresponding to the Chern field of covari 
ant vectors. The difference cocycle d(F, G) determines exactly the class 
S(a), since the coefficient of d(F, G) for a simplex a; is equal to the intersec- 
tion number of s(a@) with oj. Using these facts, Theorem I follows from 
formula (19) and Satz V. in my paper, Ref. 2. The theory on line 
elements given there holds, of course, also for covariant fields. The pas- 
sage from contra- to covariance consists merely in a factor —1 for the odd 
dimensional obstructions. 

Remarks: \. It is clear how to interpret formula (19) in Ref. 2 for the 
case of two arbitrary simple differential forms a and 8, if one observes that 
d(F, G) S(a) — S(p). 

2. Formula (20) in Ref. 2 would describe the singularities of the Jacob- 
ian of S a. 

3. If there exists a differential form of the first kind on M" then the 
corresponding line-element field carries a Chern field. See also the paper 
by Hodge,* No. 7. 

Proof of Theorem 11: We construct on 1" a special field of line-elements 
by imbedding |" in P"*! and intersecting |” by a pencil of 2-dimensional 
planes about a straight line in P"*!. For this special field one realizes 
easily that the system of ‘‘secondary cocycles,”’ Ref. 2, p. 227, consists of 
{IIj, ., I,}. Taking this field as a reference field for a Chern 
field, we get (2) and (3) again from Ref. 2, (19) and (20), respectively, ob- 
serving that the difference cocycle is II, and that one has to pass to covari- 
ance, 

Remark: Should it not be possible to imbed |" in P"*! without singulari- 





VoL. 38, 1952 MATHEMATICS: J. H. SAMPSON 895 


ties, but in a higher dimensional P” (m being the lowest dimension for 
which this is possible), then Theorem II still holds, if one interprets Il, as 
the dual cocycle class of the intersection cycle of a(k + m — n — 1)-dimen 
sional plane in P” with V". 

1 Die Idee der Riemannschen Flaeche, Weyl, H., Berlin (1913) 

? “Ueber Schnittflaechen in speziellen Faserungen und Felder reeller und komplexer 
Linienelemente,’’ Kundert, E. G., Ann. Math., 54, 215-246 (1951) 

3 “The Topological Invariants of Algebraic Varieties,’’ Hodge, W. V. D., Proceedings 
of the International Congress of Mathematics in Cambridge, 1950 


A NOTE ON AUTOMORPIIIC VARIETIES* 
By J. H. SAMPSON 
MATHEMATICS DEPARTMENT, MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
Communicated by S. Bochner, August 29, 1952 

In this note we shall make a very limited application of Theorem 1 be 
low—-namely, to the proof of the fact that compact automorphic varieties 
of the type considered admit only a finite number of analytic homeomor- 
phisms. 

Let D be a bounded domain in the space of m complex variables 2, ..., 
z,, and let fF = |v} be a discrete group of analytic homeomorphisms of D 
on itself. The elements of I’ are necessarily countable, and we shall de- 
note them by yo, ¥1, yz, ..-, Where yo is the identity element. If we write 
J (Z) for the Jacobian of y,, that is, 


O(y,Z) 


JZ) mena 
O(Z) 


then the series 


> J(Z)|? (2) 


1=0 
is majorized on every closed subset of D by some convergent series of con- 
stants.' Hence, if //(Z) is any bounded holomorphic function on D, then 


the Poincaré series 


O(Z) = So M(yZ)-[Ji(Z)}' (3) 


t=O 


also represents a holomorphic furiction on D, provided k > 2; and it has 


the periodicity property 


O(7,Z) = O(Z)-[J(Z)]—*, (4) 
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by virtue of which it is called an automorphic form of weight k. All the 
forms of a given weight & constitute a linear space, of dimension P,, say, 
which in general may be infinite. We wish to establish the preliminary 
result that in any case this dimension P, tends to infinity with k. 

To show this let A = (a, ..., a@,) be a point of D which is not a fixed 
point of any y, other than the identity yo (the set of points not fulfilling 
this requirement has complex dimension at most x — 1); and let a, a, 

, a, be any complex numbers, € some small fixed positive number. We 
propose to show that for all sufficiently large weights k there exist Poincaré 


series (3) satisfying* 


a) 
O(A) : O(A) - aj < € 
OZ; 


Denote by S the polycylinder | 2; — a; < p,j = 1, ...,n, p being taken 
small enough that S is entirely contained in D. By the remarks above we 
may find an integer g such that the quantities 


c, = sup| J;(2)| (6) 


are less than '/+ for 7 >g, the series )°c,? being, of course, convergent. For 
I1(Z) choose any fixed function, a polynomial, say, such that //(A) = a, 


- oO 
H(y,A) = Ofori = 1, ,q, and such that —-J/(A) = a, 2 H(Z) = 0 
24 02; 
at Z = 7,A for7 = 1, ,nandi = 1,...,q. Such a choice is always 
possible, since y,A # A fori = 1,2,.... Then 


O(A) R II(y,A)-[Ji(A)]}* = a + F(A), 


iaqt+l 
where 


F(Z) = > WM(yZ)-[J&Z)}*. 


i=q+l 
For Z in S we have 
“ae ee ‘ 
cf = 27°C: > (2c,)*, 
q+1 q+ 


C denoting a bound for | H/(Z)| on D. But for i > q we have 2c, < 1, so 
that (2c,)* S (2¢)*, if k > 2. It follows that, for Z in S, the inequality 


above may be replaced by 


| F(Z)| < 2-*C- 2c,)? = 2-*C (Z in S), (7) 
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the constant C’ depending only on H7(Z) and S. Obviously, we shall have 
|@(A) — a| < for large k. 


A simple computation shows that 
re) aor , 
O(A) = a +> F(A) ($ wi, ..:.9}, 


O02; 24 
and equation 7 with the iterated Cauchy formula applied to S gives us 


ma 
2 F(A)\ < ; (j = 1,..., 1). (8) 
02; 2* p a 
Thus our function 0(Z), for large k, satisfies (5). 
Then it is clear that, for large k, we can find m different Poincaré series 
O,, ..., O,, all of the same weight k, corresponding to m given systems of 
numbers (a,;, ...,; Qn), v = 1, ..., 2, such that 


i a — a,j <e€ Og = 1,:...p mp 
If € is suitably small and if det (a,;) # 0, then these n forms of weight k 
must be linearly independent. Hence, Py, > m for large k. But the 
theorem above expressed by (5) manifestly has a straightforward generali- 
zation for partial derivatives of any bounded order,*® with which a similar 
argument gives immediately 

THEOREM |. The k-genera P, tend to infinitely with k. 

The space 4 = D(mod I) can always be regarded as a complex analytic 
manifold, and we shall suppose in what follows that A/ is compact. We 
shall use the theorem just established to prove 

THEOREM 2. The total group G of analytic homeomorphisms of the auto- 
mor phic variety M on itself is always finite.‘ 

The link between Theorem | and the present problem is afforded by the 
fact that G must be a complex Lie group——possibly zero dimensional, in 
which case it is totally discrete. But G is easily seen to be compact, since 
D is bounded‘ and therefore if it 1s zero dimensional, it must be finite. 

Let di, ..., ¢», be a basis for the holomorphic densities of weight k on 
M. Since the latter is compact the P, are all finite and coincide with the 
numbers introduced above. For clearly every such density ¢ of weight k 
is the projection of an automorphic form @ on D of weight k (not neces- 
sarily a Poincaré series), and conversely. Now any transformation g of G 
transforms the ¢, into new densities $f, and we shall have relations of the 
form 


Of = Aylg)-di + ... + AuPe(g) Oe = (U = 1... Pes (9) 


the \,, depending only on g, of which they would be holomorphic functions 
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when considered as functions defined on G if the latter were positive di- 
mensional. But since G is compact, the \,,(g) would then be constants. 
We show that this is impossible for large k. 

Let ti, , |, be complex coordinates in a neighborhood of the identity 
of G, with all /; = 0 at the identity; a is of course the dimension of G, as- 
sumed positive. Differentiating equation 9 with respect to ¢; we have 


O 


d;, 
a,” 


and at 4; , - = () this is 


n 
Op o 4 
D> — "05 = 0, (10) 
a=] 02, 
45 Oz, | es Dos , 
where the 7,4) are, of course, the infinitesimal generators of G, 
ty 0 
, 2, being local coordinates on M. Since the 5) cannot vanish 


lo 
“ls 


everywhere on AM/, we can by Theorem | choose k& sufficiently high that 


equation 10 cannot hold. 


* Written while the author was under contract to the Office of Naval Research. 

! Siegel, C. L., Analytic Functions of Several Complex Variables, Inst. for Advanced 
Study, Princeton, 1949, pp. 124f 

2 Giraud, G., Legons sur les fonctions automorphes, Gauthier-Villars, Paris, 1926, pp 
19-25 

§ Giraud, G., loc cit. 

4 Andreotti, Aldo, Sopra il problema dell’uniformizzazione per alcune classi di super 
ficie algebriche, Rend. accad. Naz. XL, Series IV, I, 16 (1951) 

Bochner, S., “On Compact Complex Manifolds,”’ J. Ind. Math. Soc., XI, Nos. 1 and 2, 
1-21 (1947). 

Hawley, N., “A Theorem on Compact Complex Manifolds,” Ann. Math., 52, No. 3, 
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BARNARD’'S CSM TEST OF SIGNIFICANCE 
By EpwIn B. WILSON 
OFFICE OF NAVAL RESEARCH, BOSTON 


Communicated August 21, 1952 


As Pearson remarked! in following one article by Barnard? and preceding 
two others:* “The problem of testing the significance of a difference be- 
tween two proportions is one which receives early attention in textbooks 
on mathematical statistics, and it might be thought to be one of the ques- 
tions whose final solution lies behind us."’ It is my impression that in many 
textbooks the problem receives very little attention; a rule is stated with 
small discussion of its possible “‘outs’’ to contrast with its ‘‘ins,’’ and the 
student is indoctrinated into its mechanical use by example and exercise. 
If there were still basic aspects of the problem to discuss in such a journal as 
Biometrika only five years ago, how long will it be before the flood of new 
textbooks which constantly inundates the market will be giving an ade- 
quate account of the problem ? 

Barnard” performed a very real service to elementary statistics in dis 
cussing different sampling procedures which might lead to a fourfold 
table and some of the different tests of significance appropriate to them. 
Barnard® * postulates that the significance of a result is to be judged rela- 
tive to the set of all possible results which might have been obtained by the 
same procedure applied under the same external conditions of the sampling 
process employed, and in particular that the significance level P of a par- 
ticular result is the sum of the probability of that result and of all those 


equally or more unfavorable. It remains to define in a particular scientific 


problem, or for certain types of problems, what is to be meant by ‘“‘equally 
or more unfavorable,” and as there is still difference of opinion as to proper 
fundamental postulates, one can only expect that there will be difference 
with respect to this subsidiary definition. 

How difficult almost to impossibility it is to state one’s postulates and 
definitions clearly, and without possibility of misleading oneself or one’s 
readers with connotations which go beyond that which one would denote, 
is familiar to all students of the foundations of mathematics, and even more 
to many others struggling with less definite fields of endeavor involving 
communication with one’s self or others. This may be illustrated by refer 
ence to my initial quotation from Pearson in which he writes of the problem 


‘ 


of testing the “‘significance of a difference between two proportions.” If 
we have from one universe a sample of m with a “successes,” and from 
another a sample of m with c, the observed difference of the proportions is 
Pa — Pe = a/m — c/n, Inthe set of (m + 1) times (nm + 1) results which 


might arise, the (m + 1) times (7 + 1) differences p, — p, are arranged by 
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their numerical values in a serial order from —1 to 1, and if we are con- 
cerned only with the significance of the difference itself, we are over one 
very hard hurdle in that we have equal or greater numerical values of the 
difference to define “equally or more unfavorable.” 

Now Pearson may mean this and nothing more, but I believe he is of the 
opinion that equal differences are not of equal significance in judging the 
significance of a particular result obtained in the comparative experiment 
so frequently used in science. This may be quibbling, but I think not; 
and it certainly is not intended as any criticism of Pearson, but merely as 
an illustration of a difficulty of perfectly precise expression. There is the 
obverse difficulty that, in the effort for precision, expression may become 
so involved as to become unclear for that reason! The greatest importance 
of the distinction drawn in my excursion into semantics is, however, the 
emphasis I hope thereby to throw on Barnard’s achievement of a novel 
method of arranging the (m + 1)(m + 1) results which can arise in a com- 
parative experiment in a serial order without the use of the difference 
Pa — Pp, or of p, or p, at all, and thereby getting over the first very difficult 
hurdle toward the determination of the significance of the result observed 
in such an experiment. 


The observations may be represented by the paira,b; c,d, where a + b = 
mandc +d =n. It is considered that the two pairs a,b; c,d and 6,a; 
d,c are equally significant with respect to the “null hypothesis’ that the 


pairs in the set arise from a pair of universes in which the probability p is 
the same and the result a,b; c,d or its symmetrical b,a; d,c had arisen 
merely by the chances of sampling. Then a start is made at the two ex- 
tremes m,0; O,n and O.m; 1,0, to which taken together the significance 
level P = p"g" + q"p" is assigned.’ Next one considers the two adjacent 
results m — 1,1; Om and m,0; 1,n — 1, and of course their symmetricals; 
one writes down the two alternative significance levels 


P = p™q" + qg™p" + m(p™—'q*t! + g™—'p*t!) and 
P’ . p™q" + q™p” + n(p™ tla" 1 +. qut'p™ 1) 


and one determines as next in order of significance that one of the two 
which has the smaller maximum, and so on.® For the details reference 
must be made to the origifal, but at least this account of the method 
seemed necessary to make at all clear the points I would make. 

Barnard gives the table form = n = 7 arranged in the order determined 
by his method. It is obvious that his serial order is different from that 
based on the difference p, — p,. He gives the values of the significance 
levels P at maximum according to his method and also the values which he 
assigns to the same results by R. A. Fisher’s so-called exact method,® and 
it is seen that his serial order and Fisher’s differ, and that his maximum 
values are considerably lower than Fisher's, and indeed in most of the 





VoL. 38, 1952 MATHEMATICS: E. B. WILSON 901 


cases I have tried, Barnard’s method gives significant estimates very much 
closer to the old-fashioned estimates based on pz — p,. If time should 
validate Barnard’s method or any which gave resulting significances ap- 
proximating his, it would also support the opinion I have long held that 
Fisher’s ‘“‘exact method”’ is in fact very inexact and very wasteful of sig- 
nificance in small samples—-a denouement most ironical in view of the em- 
phasis Fisher has so properly and effectively placed upon the desirability of 
extracting from data all the information which they have in them, 

Many years ago when Yule became suspicious of the rule for estimating 
the degrees of freedom in the application of x*, he succeeded in setting up a 
game of chance whereby he showed, at least in that case, the incorrectness 
of the rule under suspicion and could discern the correct rule.’ I have in 
the past dozen years given much thought to designing a simple chance ex- 
periment to do something similarly effective for this problem of significance 
in small comparative experiments, but without success. As early as 1936, 
the late Barbara S. Burks told me that from her work with large numbers 
of very small samples of twins and of identical twins she was convinced 
that the exact method seriously underestimated the significance level. 
She was a good statistician fully aware of the dangerous bias, the subtle 
misjudgment which can so readily deceive one in such an inference. I 


hoped she would print her findings on this methodological point, but I be- 


lieve she died without doing so. 

There is another remark one may make on the tables given by Barnard, 
namely, that the (maximum) significance of two tables with rows and col- 
umns interchanged need not be the same. Thus 7,2; 0,5 and 7,0; 2,5, 
the former obtained by sampling one universe by drawing 9 and the other 
universe by drawing 5, the latter by drawing 7 from each of two universes, 
are given the respective significances P = 0.0062 and P = 0.0070. This 
is so interesting an observation that I wish to illustrate it on two simple 
cases. 

Case 1. Consider the array of eight tables for a sample of 3 and a sample of 1, and 
also the array of nine tables for a pair of samples of two each, namely 


3,0 


The result 2,1; 0,1 occurs in the former and the result 2,0; 1,1 with rows and columns 
interchanged in the latter. It is easily seen by the simplicity of the arrays that 


P = pq + 4p2q? = x + 4x?, P’ = 4p’ '— 6p'%q'? = 4x’ — 6x’? 
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are the two respective significance functions of the two tables, i.e., they give as a function 
of the unknown hypothetical proportions p or p’ common to the respective pairs of unt- 
verses the total probability of the given tables and of all equally or more significant ones 
in its own array. Note that the maximum of P is 0.5 at p = 0.5, which is the marginal 
value seen in the table, and that the maximum of P?’ is 0.625 at p’ = 0.5 also, as it hap- 
pens, which is not at all the value (0.75) which would be read from the table. 

It follows that for all values of p’ which make P’ larger than 0.5, i.e., for p’g’ > 1/6 
or for p’ between about 0.788 and about 0.212, there can be no value of p for which ? = 
P’. {tis clear that the variables most simply involved are x = pq and x’ = p’q’, either 
of which may in this simple case be expressed in terms of the other if we impose the 


equality 7? P’, namely as: 
Y= /A(l- V1 x — 6x2), : + V/1 + 64x’ — 96x’), x’ S 1/5 


For each value of x or x’ there will be two values of p or p’. The following table may be 
worked out 
0O9or 0.1 O8Sor0.2. 0 7o0r0.8 0. 60r0.4 0.5 
0.967 0.919 0. 864 0.8138 0.788 
0.033 0.081 0.136 0.187 0.212 
0.1224 Q. 2624 0. 3864 0.4704 0. 5000 
Cast 2. If we take the results 7,1; 0,6 and 7,0; 1,6, for which Barnard’s tables (p. 
137) give P? Pp’ 0.0018, being two significant figures for 15/8192, which is attained 
for p = p’ = 0.5, so that the maximum values of ? and P’ coincide in position and in 
amount, it is still true that the functions ? and P’ are not identical.2 As a consequence 
it is still true that values of x and x’ must occur in pairs if the condition ? = P’ is im 
posed, and to each x and its correlative x’ there will be a pair of values of p with their 


correlative pair of values of p’; the equations are: 


P 6x 23x* + 26x7 and P’ 14x’6 26x". 


It is clear enough from these cases that the interchange of rows and 
columns in a fourfold table, which makes no difference in its significance 
by Fisher's exact method or by x? (or by testing py — p, against its standard 
deviation as formed from the margins, which is equivalent to x*), does 
make a difference with Barnard’s method unless one is willing to assign to 
the unknown common probabilities some one of an infinite number of 


pairs of values determined by the equation P = P’. The more I reflect 


on a considerable number of such cases, the less I like this idea. I should 
perhaps prefer to admit that the significance of two tables which differed 
only by interchange of rows and columns might thereby differ. But I 
find in the literature :’ 

“The marginal totals. ..represent two independent dichotomies of the 
data; they do not constitute a double dichotomy. The problem may be 
put in either of two ways: (1) Is the division between the effects the same 
for both methods? (2) Is the division between the methods the same for 
both effects? These questions are logically identical, though for a particu- 
lar problem one of them seems more natural than the other.”’ 

I should readily admit the truth of this statement if the division were a 
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double dichotomy, but have never been sure of it for the comparative ex- 
periment. If it is true, I suppose Barnard’s method will have to go over- 
board. 

One of the best methods of thinking over a problem when no particular 
solution thereof is commonly accepted is to work through a variety of 


special cases. It may save some time of those who may still be wonder 


ing, as | am, about this matter of the comparison of two percentages or 
proportions if I should give the details I have found for two of the sim- 
plest of the several cases which are contained in my dossier on this sub 
ject. I shall choose the case of two samples of three and that of one 
sample of four and one of three. 


Two Samples of Three.—-One starts with the result 3,0; 0,8 and its symmetrical for 
which together the probability is 2p%q° which cannot exceed !/j. or 0.031 no matter what 
the unknown common probability assumed and thus fixes Barnard’s significance as 
0.031. One goes on to the obviously equally significant pairs 2,1; 0,8 and 38,0; 1,2 
taken together, with their symmetricals be it understood, and have a total cumulated 
probability 6x? — 10x° with maximum ' =(.219. (One may check that this maximum 
is less than that, °/16, which would have been found if one had tried to put 2,1; 1,2 in 
second place.) In the third place one has to try the pairs 1,2; 0,5 and 3,0; 2,1 against 
2,1; 1,2. The respective probability functions are 6x 10x? + 2x? and 6x? + &x3 
The respective maxima are 0.5097 and 0.5000. — It is a close decision but 2,1; 1,2 wins.*? 
Then we have for 1,2; 0,3 in fourth place P = 6x 18x? + 20x? with maximum 0.6875 
Next comes 1,2; 1,2, one of four equivalents, with cumulative probability 6x Vx? + 
2x8 and maximum 0.969 

We may summarize these findings and some that may offer additional food for thought 
in the following table. The columns are: (1) serial order, (2) ‘‘result,’’ (3) funetion P, 
(4) its maximum, (5) x therefor, (6) value of x from table, (7) pa p-, (8) maximum P 
for serial order determined by Pa pe, (9) Barnard’s value for Fisher's exact method as 
I interpreted it, (10) Fisher's one-sided value as insisted on by him in Serence, loc. cit.® 


(3) 4 5 { ) bad ) 10) 
0.031 0.25 0.031 0.05 
0.219 0.25 ) G.219 0.20 
0.500 0 ) ; 0 OSS 0.50 
0 688 0.25 0 688 0 50 
0.969 0.25 22 ? 
I l 


In comment I may say that I do not like Barnard’s method in respect to its giving 
0.969 instead of 1 for the significance level in V(4) because it seems to me that as the 
result 1,2; 1,2 shows no difference in the observed proportions, the value should be 
strictly 1. Further, the values attributed to the exact method by Barnard (as inter 
preted by me) in column (9) show the unsoundness of that method as thus interpreted 
And in column (10) I have had to put ‘‘?”’ opposite V because I do not know how Fisher 
would sum up one side of the series 0.2, 0.6, 0.2—perhaps he would get 0.8. It is note 
worthy that as shown in (5) the maximum occurs for all results of the comparative ex 
periment with two samples of three at the value x = 0.25, i.e., for p = '/, although for 
result IV the tabular value is p = '/s, Le., x = O14. 

One Sample of Four and One of Three Phe findings for this case may be tabulated 
with the same conventions as before 
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(2) 7 (4 5 +f (8) (9) (10) 

4,0;,0,3 x 0.016 016 029 0.029 

4D: 1,3 3x? 0.062 125 143 0.143 

3,1;0,3 3x? 0.125 078 143 0.113 

0.218 153 429 0.429 

0.305 411 429 0.286 

0.453 312 486 0.371 

0.563 563 1 0.571 

0.797 797 1 0.629 

1X 3,1; 3 3 0.G84 O84 1 0.764 
X -§ l 0 1 1 1 


to to = tt te 
- UN oe BS oS 


By looking at (1) and (7) one perceives the very different serial orders set up by Bar- 
nard’s method and the old method based on p» p:,. By comparing (4) and (8) one 
sees the very different significance levels at maximum which are assigned to the same 
tables by the two methods. Chi-square would give a still different order, viz., I, III, II, 
V,1V, VI, VII, VIII, IX, X; corrections can be applied to chi-square but one cannot 
expect them to be very good in such a case. Indeed until we know more definitely what 
is the proper value to assign (not merely one of several possible maximum values), we 
have no way of telling how good any such correction is. 


1 Pearson, E. S., Biometrika, 34, 139-167 (1947) 

2 Barnard, G. A., /bid., 123-138 

§ Barnard, G. A., [bid., (a) 168-169; (b) 179-182 

‘Ifm+n2=(m n)*, there will be a single maximum at p = g = '/, of amount 
Pawl/ymtr-l, lfm+n<(m n)*, i.e., if m being assumed the greater, we have m > 
n-+'!/,+! VW 8n +1, then there will be two maxima in this extreme case and it is 
not easy to solve for the values of p org = 1 p 

5 Owing tothe symmetry of ? if pis a value giving a maximum, so is q; the expression 
for P may be converted into one in the new variable x = pq. This cuts in half the order 
of the algebraic equations with which one must work. The values of p and g may be ob- 
tained from that of x by 


p lt V1 4x), q = '/(1 F V1 — 4x) 


Moreover the probability of a pair of symmetrical results always appears as the product 


s : k k . ‘ 
of a power of x by the sum S, = p* + q°, and this sum can immediately be expressed as 
a polynomial in 1 tx, with no radicals, There is, however, a simple recursion formula 


Se+1 = S, v.S, 1, Which permits one to run up a set of .S’s very easily, starting off as: 
So, = 2,5, = 1,5. = 1 2x, Sy; = 1 3x, I have found these devices useful in con 
ducting a considerable long-hand exploration of Barnard’s method; but there may be 
even better devices 

6 Barnard says that these are the significance levels Fisher would get. He further 
states (p. 124) ‘The probability (1) added to those of all results less probable than that 
is the level of significance according to Fisher's exact method. It was my impres- 
sion that one had to add not only all the probabilities which are less but also the one 
that is equal if there is one, and so far as I have checked his values that seems to be 
what he actually has done. (I believe the four values 0.283 which are entered near the 
N.W. and S.E. corners of the table were accidentally not doubled and should read 0.47.) 
But I cannot be sure that Fisher would agree with him. For I applied identically that 
rule in Science, 93, 557-560 (1941) toa result carefully stated so as to be purely compara- 
tive and completely symmetrical, and Fisher promptly recomputed the example, Jbid., 
94, 210-211 (1941), adding up only one side of the distribution and thus getting only 
one-half of my values. As I had also chosen my example so that adding up both sides 
would give more than 6.05 but adding up one side alone would give less than 0.05, I came 


aR ga Sa 
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to the conclusion that the result was not significant, whereas Fisher pronounced it sig- 
nificant. Apparently either Fisher changed his mind as to his exact method between 
1941 and 1947, or Barnard is wrong in his statement of the method, or perchance Fisher 
was discussing a different problem from the one I set. 

7 Yule, G. U., J. Roy. Statist. Soc., 85, 95 (1922) 

8 The simplest case of this sort that I have found is 3,0; 1,2 and 3,1; 0,2 

® Eisenhart, C., Hastay, M. W., and Wallis, W. A., Selected Techniques of Statistical 
Analysis for Scientific and Industrial Research and Production Management Engineering, 
McGraw-Hill, 1947, Chapter VII. 

10 In the early part of his paper Barnard says that what one should try to do is to make 
the probability ? as nearly constant as he may throughout the interval p = Oto p = 1, 
i.e., to make P as independent of p as possible. Then he actually chooses the lower 
maximum as the criterion for ordering. It is obvious that although this criterion may 
in most cases actually select the probability function which is most constant. it has 
nothing to do logically with any which might be based on the effort to get constancy. As 
Barnard clearly wishes to play safe with respect to errors of the first kind at the expense 
of errors of the second kind (since he determines maxima of P), I should suppose that 
one way to formulate the condition of maximum constancy would be to use the difference 


+] ; : . ; ; 
D = Pu. — - P dp or its value D/P yx relative to the maximum. I raise this sub- 


ject at this point because if we were to examine the deviation D of the function 6x 
10x? + 2x4 for 1,2; 0,3 from its maximum, we should find 0.0954; whereas for the fune- 


tion 6x? + Sx for 2,1; 1,2 the deviation from its maximum is the much larger value 
0.2429. Thus upon this criterion the decision would not be close and would place 1,2; 
0,3 third with the significance level 0.510, leaving 2,1; 1,2 to come in fourth at the level 
0.688. This is the kind of instability ene must expect when dealing with small numbers 
and illustrates how critically significance levels in these cases depend on our assump- 
tions and criteria and how careful we must be in stating them if we would really discuss 
the problem of comparing two percentages, 


THE NEWTONIAN POTENTIAL OF A SPHERE OF A EUCLIDEAN 
SPACE Fy EMBEDDED IN A EUCLIDEAN UNIVERSE E,* 


By EDWARD KASNER AND JOHN De Cicco 


DEPARTMENTS OF MATHEMATICS, COLUMBIA UNIVERSITY, NEW YORK AND DEPAUL 
UNIVERSITY, CHICAGO 


Communicated August 2, 1952 


1. Consider the two forces of attraction fF, and Fy exerted by a homo- 
geneous material (V — 1) dimensional surface SS and NV dimensional volume 
V of a sphere of radius r and mass V7 of a Euclidean space /}y upon a 
particle P of mass m in /y at a distance a from the center of the sphere 
such that /y is embedded in a Euclidean universe /, of n dimensions where 
2< N <n. Both of these two forces are directed along the line joining 
the position of the particle P to the center of the sphere. These two fields 
of force are conservative and solenoidal.! 
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11 Td 


: sin’ 2 0(1 — b cos 0) dé 
I = I(N, n, 6) ‘ (1) 
0 (1+ 4? — 2b cos 6)"/? 


where b = r/a, then the force of attraction F., exerted by the material surface 
S of the sphere on the particle P is 


kmM I(N, n, b) 


rt N-—- 1 
B ( 
o “2 


and the force of attraction Fy exerted by the material volume V of the sphere 


on the particle P 1s 


kmM N : ; p 
Fy, : : pv tN, n, dp. (33) 
Ite N DTA a 
PB “eae! 


These expressions are deduced by means of spherical coordinates in /2y 
from an extended form of Newton’s Law of Universal Gravitation which 
states that the force of attraction between two particles in a Euclidean 
universe /¢, of mn dimensions 1s directed along the line joining the two par- 
ticles and varies directly as the product of their masses, and inversely as 
the (n |) power of the distance between them. 

1.2. If the function J is detined by the integral 


n—3 


: pe A 
J(N, n, X) 


(A + @) 2 


where \ = (1 — b)*/4b = O, so that J obeys the condition 
V-n-1 


J(N, n + 2, X) IA(1 + A)] 2 J(N, 2N — n, d), 


then the integral I is 


, bICN, nd) + (1 — b2) ICN, n + 2, d)). 


» 


N+3 h* 


(0) 


If in the preceding formulae, we substitute V = n, the following results 


(1.3) and (1.4) are obtained. 

1.3. The force of attraction exerted by a homogeneous material (n — 1) 
dimensional surface S of a sphere of radius r and mass M upon a particle P 
of mass m ata distance a from the center of the sphere is directed along the 
line joining the position of P and the center of the sphere. Its value is 
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kmM |. 
u 


n 1 


FPF, = kml 


1 
,% 


2a" 


0 if 


When the particle P is in the exterior region of the sphere, the force 1s 
such that the mass ./ of the surface S of the sphere is concentrated at its 
center, and when the particle P is in the interior region of the sphere, the 
force is zero. 

1.4. The force of attraction exerted by a homogeneous material n dimen 
stonal volume V of a sphere of radius rand mass M upon a particle P of mass 
m at a distance a from the center of the sphere is directed along the line joining 
the position of P and the center of the sphere. Its value ts 


1} 


| km\M 


a 


| kmM 


n 


F, 
l 
r 
When the particle P is in the exterior region of the sphere, the force is 
such that the mass \/ of the volume V of the sphere is concentrated at its 
center, and when the particle P is on the surface or in the interior region of 
the sphere, the force 1s proportional to tne distance a of P from the center 
of the sphere.” 
2. The integral representation of the hypergeometric function is 


I'(y) “1 
F(a, B, ¥, 3) ere / -'1 — t)y-F- "1 — at)" dt, (9) 
L(B)V(y Dp) 0 


where R(y) > R(8) > 0. This is uniformly convergent in any closed do 
main of the complex z-plane provided that the z-plane is cut along the real 
axis from +1 to ©. This analytic function has branch points ats = +1 
and zg = o,? 


If m and N are positive integers such that 2 < N n, set 


;, | eet ee 
Gla, 2) hk (a, a 4 N i £e 


By (4), (9), and (10), it 1s found that 


(“ ye ) 
B t 
2 2 


[n 
J(N, n, X) —— G = 
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Upon substituting (11) into (6), the result is 


: ; y F a 1b 
I(N, n, 6) - 4 
211 2 2 (1 — bd)? 


| - bd) 2N—n-—2 1b 
_- » (12) 
( 1+5 2 (l — >) 


where the plus or minus sign is taken according as0 Sb < lorl<b< 
4 o, 

2.1. In terms of hypergeometric functions, the force of attraction Fs of 
(2) 1s given as follows. If0<r<a< +, then 


: kmM ; v4 lar 
Fs au ,- ,) + 
(a — r)? 


2a(a — + ; 2 


N 2N — —2 | 
 ( <S  — )}. (13) 
) (a — r)? 
+o, then 


kmM ja (" ; lar ) 
a)"?L~ (a — r)* 


(: - *)" C (2 —-n-—2 4ar )| (14 
r+a : 42 : (a — r)? ; ” 


From (7) and (13), it is deduced that 


har 
l, - ~) = 
(a — r)* 


r n-—2 
(1 _ ) if0<r<a<+eo,. (15) 


a 


2.2. In terms of hypergeometric functions, the force of attraction Fy of 
(3) ws given as follows. IfO<r<a< +, then 
£ 


; kmMN [" p* : , lap 
fk \ i ‘ / ” ’ pai .) + 
ZarN Jo (a—- op : . (a — p)* 


a - eee RN. - § fap 
G -}|dp. (16) 
atop y (a — p)* 


lffosa<r<+ o, then 


kmMN [*" 
diel OnyN 
Zar 0 (a 
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lap kmMN 
~} | 20 — x 
— p)* 2ar% 
J . m lap ) 
a+ (p — a)"~? y (a — p)* 
p—a\"* /(2N - — : lap 
G a dp. 
pt+a : (a — p)’ 


From (S) and (16), it is deduced that 


r n—l ‘ 
p Mk, = hap 
f , ( ) - 
0 (a — p)" 7 ; = p)* 


Also by (8) and (17), it is found that 


*a- n—] ‘ 
p tn — lap 
0 (a: — p)*"* y y (a — p)’ 


n#O0S a< 
n 
This is also a consequence of (18) by letting r approach a. 

3. Let U(x) denote the Newtonian potential of the field of force pro- 
duced by the attraction of a body Vy upon a particle P of mass m in a 
Euclidean universe /, of m dimensions where | S N Sn. The rectangu- 
lar components of the force vector are X; = ax! for? = 1, ..., mM. 

x 
3.1. The Newtonian potential U(x) of the field of force Fs of (7) is given 
as follows. If n 2 3, then 
kmM 
(n — 2) a” 
l 
kmM 
(n i 


ee) Or: 


If-n = 2, then 


l 
Jamar log 
a 
l l — 


] 
[ema log oY: 0:3: 6-8°8: 
r 


For the potential function l’ of (20) and (21), it is found that V°?U’ = 0 
provided that a # rwhereO0<r<+0,0S5a<+o., 
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3.2, The Newtonian potential L(x) of the field of force Fy of (&) is given 
as follows. Ifn 2 3, then 


SE ee <a 
~ 2) a” 


(22) 


km M a? 2 az 
| ( + Sasr<+o, 


n ») 


r 


kmM 
| (n 


Ifn 2, then 


| 
Jemar log if tri<a< +o, 
( 


l 
kmM a” yr? 1 ~ 
+ + r® log Sasr<+4+o. 
ea 2 2 r 


The Laplacean of the potential function Ll’ of (22) and (23) is 


| 0) au O<r< ak ao, 


| a 
Vil — kmp . if OS aSr0<Kr< te, (24) 
) 


| 


where yw is the constant density. 

1. In this final section, the potential function (’ of a homogeneous ma 
terial surface of a sphere of an N dimensional Euclidean space embedded 
in a Euclidean universe /, of m dimensions where 2 <= N S nand3 Sn, 
is discussed. ‘This is the Newtonian potential l’ of the field of force Fs. 
In the special case where the particle P is in /y, the force Fs is given by 
(2). 

4.1. If pand qare the minimum and maximum distances of the particle P 
from the surface S of the sphere, then 


km AM | " sin’ *@d0 
N-1 1 A 4 e\2—= (25 
(n 2) B( , f ( : ) 2 (<9) 
2 2 


2 O¢ s? g: SI 2 
p* co _ os 111 2 


The particle P is in the Euclidean space which is completely orthogonal 
to /-y if and only if it is at a constant distance a from the surface S of the 
sphere, in which case the potential is LU’ kmM/(n — 2)a"*. As the 
particle P approaches the surface S of the sphere, the potential (’ remains 
finite if and only if either V n or else both N and n are even integers 
such that 2 $ N <n.‘ 

Pe St ee ' ; 
b2.. dt 2:3 - <= N <n, let U, and U2 denote the potential functions 


(25) relative to the euclidean universes I.,, and Ey y42, respectively. Then 
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2N—n ; 
- x : (26) 
2 (mn ——.2)¢be)e" 
This is deduced from the fact that, by (4) and (5), the potential Ll’ of 
(25) can be written as 


kmM . UN ' 
2° JUV, 2, A) 
Se ‘. 


n - 


(n — 2)(q* — p*) * 


n-2N 


_ kmM(q? — p*) ? 7 
(n — 2)(pq)"—*— a( 
where \ = p?/(q? — p*) > 0. 


4.3. If 0 <p <q, then in terms of hypergeometric functions, the potential 
function (25) ts 


kmM G(" -2 r) 
(n — 2)p*-* : a p* 


kmM 


(n — 2)p* 


From (28) is deduced the formula 


nN /2N N-1  — g? 
(’) i ry ht ) (29) 
q 2 2 p* 


If the particle P varies over the (n — N) dimensional surface S; of a 


sphere in a Euclidean space /:,.y 4; such that S») 1s completely orthogonal to 
the Euclidean space /y which contains the surface S, of the original sphere 
in y, and S, intersects /y in two points inverse with respect to S;, then 
the potential U of P is inversely proportional to the (n — 2) power of the 
minimum distance p (or maximum distance q) of P from 5). 


* Presented to the American Mathematical Society, September 1952 

1 Kasner and De Cicco, ‘Potential Theory in Space of m Dimensions (Part 1),’’ Proc 
Nati. ACAD. Scr., 38, 145-148 (1952) 

2 The formulae (7) and (8) are extensions to m dimensions of the corresponding classical 
three-dimensional cases. See O. D. Kellogg, Foundations of Potential Theory, Frederick 
Unger Publishing Co., New York, 1929 

3 Whittaker and Watson, Modern Analysis, Cambridge University Press, 1950 

4 This result was obtained by Tom Arai and Charles Stewart, graduate students at De 
Paul University. 
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THE COPPER METABOLISM OF DROSOPHILA* 


By D. F, Poutson, V. T. BOwEN, R. M. HILsE AND A. C. RUBINSON 
YALE UNIVERSITY AND BROOKHAVEN NATIONAL LABORATORY 


Read before the Academy November 2, 1951; received July 29, 1952 


The copper metabolisin of insects has been discussed in the past by 
various authors, mainly on the basis of analyses of whole animals. In the 
fly Lucilia cuprina, Waterhouse! studied by sensitive histochemical methods 
the uptake of copper and its distribution to the tissues. The histochemical 
findings were supplemented and supported by careful chemical analytical 
methods. 

One of us has made, and in part reported elsewhere?~* similar histo- 
chemical studies on a number of species of Drosophila. However, the work 
of Okamoto et a/.,° on copper in tissues and our knowledge of such copper 


proteins as hemocyanin® and ascorbie acid oxidase’ indicate the presence 


in tissues of copper compounds which will not react in the various histo- 
chemical tests. It thus seemed desirable to extend our observations by 
tracing the copper metabolism of Drosophila larvae with the radioisotope, 
Cu. This has the advantage of providing both a somewhat more sensi- 
tive test for copper and one which is independent of the state of combination 
of the copper in the tissues. 

Methods. The radioisotope Cu® was produced® by neutron irradiation. 
The c. Pp. metal was cleaned in advance and irradiated wrapped in cleaned 
aluminum foil. After short (36-48 hrs.) periods of irradiation, careful 
determination of decay curves of various samples has confirmed that no 
significant contamination by other radioisotopes has been produced.  Ac- 
cordingly no chemical separation of the irradiated copper was made. 

The irradiated copper, carefully weighed before irradiation, was dis- 
solved in HNOs-H.SO, and the HNO, driven off by heating to SO; fumes. 
The solution was made up to contain a known weight of copper per gram, 
and suitable aliquots were added to Drosophila medium’ or to yeast media. 
By mixing the Cu® solution into 50-100 g. of hot fly food, high uniformity 
of distribution was secured as measured either radiochemically or by chem- 
ical analysis. 

Detection of radioactivity has been by means of conventional Geiger- 
Miiller counting of whole or dissected larvae; of contact auto-radiography 
of dissections of larvae; and of histological autoradiography of paraffin 
sections. Data obtained by the two latter methods will be presented 
here only in so far as it relates to (confirms or extends) the counting obser- 
vations, and will be reported in full elsewhere. 

By parallel counts on aliquots of the original Cu® solution, and on weighed 
portions of the radioactive food, and by giving due attention to decay cor- 
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rections and to uniformities of counter response, it has been possible to 
convert much of the data from counts per minute directly to micrograms of 
the original copper sample. Although the energies of Cu® 8 particles are 
low (0.57 and 0.66 Mev max.), the very small mass of the tissues being 
counted has permitted us to neglect self-absorption corrections, and to 
obtain results of adequate precision. 

The larvae used were raised from eggs collected over 12-hour periods. 
Early third instar larvae selected for uniformity of size were used in all the 
experiments. Except during the time required for manipulation and 
transfer the larvae were kept in an incubator maintained at 25°+0.1C. 

Larvae were scoured of radioactive food by transfer for periods of 2 hours 
or more to non-radioactive medium. Before preparation for counting or 


TABLE 1 
UprakKE OF Cu BY DrosopHILA LARVAE. Errectr OF CANDIDA ALBICANS (TEXAS Y 12) 


cu®t LEVEI C./MIN 
SPECIES 3./G TIME OFF LARVA 


Repleta 2.8 § 4'/, 108 
; 3 113 

5Cu 101 

4 { 10 

38 11 

Melanogaster 2 8 3 28 
. ; 35 


Ananassae 


Notre: Similar data obtained for Saccharomyces (baker's yeast ) 


fixation larvae were thoroughly washed in Drosophila Ringer solution. 
Larvae used for whole counts and sectioning were fixed in hot formalin- 
alcohol. For contact autography larvae were dissected without fixation 
on dry slides on which the organs were spread and subsequently 
oven-dried at 60° for 30 minutes. For counts of dissected larvae the dis 
section was performed in the dishes in which counting was to be carried out 


and parts transferred to other dishes so that there was no loss of material or 


activity. 

Experiments and Results.Since a major fraction of the food of Droso- 
phila larvae consists of the yeast cells growing on the medium, it was 
first desirable to examine the effect of yeast on uptake by larvae of copper 
from the medium. The data from one experiment designed to test 
this are set out in table 1. The copper was mixed into the hot fly 
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food, in the concentrations indicated and replicate samples were taken. 
One of each pair was inoculated with yeast which was then allowed to grow 
for 12 hours. The other was left sterile. Larvae were placed on the two 
types of media for S hours and removed as indicated in the table. The uni- 
formity of counting rate of larvae taken from the same copper concentra- 
tions indicates that even though there had been very appreciable yeast 
growth on the inoculated cultures, the yeast cells have no effect on the 
amount of copper absorbed by the larvae. Similar experiments over other 
and longer times and using Saccharomyces as well as Candida, confirmed 
this. 

The effect of varying concentrations of copper in the medium on the 
amount taken up by the larvae at constant time, is illustrated in table 2. 
Although uptake ceases to be strictly proportional to concentration for 
ananassae between 10 and 25ug: Cu per gram of medium, and for repleta 
between 2.5 and 5ug. Cu per grain, the control of uptake is not strict. In 


TABLE 2 


UpraKe OF Cu® By DrosopHILA LARVAE. Errecr or Torat Cu IN THE MEDIUM 


CU CONTENT AFTER 8 HRS. ON AND 4-6 HRS. OFF, DIVIDED BY THAT 


POTAL Ct AT 2.5 wG, Cvu/G 


UG./G REPLETA MELANOGASTER ANANASSAE 
0.25 0.1 0.1 
2.5 | 
5.0 
10.0 
25.0 
100 
200 


repleta even at 200 wg. Cu per gram, there is no evidence of a real plateau 
in the curve relating uptake to concentration, With the smaller ananassae 
and melanogaster, increase has ceased and an actual decrease in uptake is 
observed between 100 wg. Cu per gram and 200 ug. per gram. Typical 
symptoms of copper toxicity make themselves evident at higher concentra- 
tions. 

A similar difference in capacity for copper uptake is seen in table 3 which 
illustrates the relation between copper content and time on the medium for 
repleta and ananassae. In these experiments larvae were allowed to 
feed on radioactive food for periods of 4 hours, 8 hours, 16 hours and 24 


hours respectively, and then on non-radioactive food for 2 hours before 
preparation for counting. It is clear that uptake at the several con- 
centrations reaches a plateau after about 16 hours on the Cu®! medium in 
both species. Accordingly 16 hours has been used as the feeding time in 


the experiments on excretion, 
Excretion of copper has been followed by determination of activities at 
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intervals after the standard time of 16 hours on food containing Cu™ at con 
centrations ranging from 2.5 to 25.0 ug. Cu/g. <A typical experiment on 
the excretion of absorbed copper is illustrated in table 4. Here two hours 
has been taken as the standard flushing out time, since it is established that 
the solid contents of the gut are wholly passed through in somewhat less 
than 2 hours. In the case of copper we do not seem to have the problem of 
unabsorbed material lying for long periods between the peritrophic mem 
brane and the gut epithelium as has been described for barium". Although 
in this experiment the ananassae gave very variable results, it is plain that 
both species do excrete copper absorbed from the medium. In repleta 
excretion proceeds more gradually and steadily to a lower level of retention 


rABLE 3 
UpraKkE oF Cu BY DROSOPHILA LARVAE 


C /MIN./LARVA AFTER T ON AND 2 OFF 
RPFECT OF TIME ON FOOD 
C./MIN. ‘LARVA AFTER 4 HRS. ON, 2 OFF 


KEPLETA ANANASSAE 


1 HRS S HRS 1} HRS 1 HAS 8 HRS 16 uRs 24 Wes 


FABLE 4 


Tr ore 


SS OF CUS* ON PLAIN FOOD Cc MIN LARVA 16 HRS. ON Crt RATIO 
AT 2 URS OFF 


POUTAL Ct REPLETA ANANASSAE 
4 HRS 8 HRS 16 HRS 1 HRS S HRS 


68 5! 0 25 0 65 0 57 
67 BY 0.93 
29 0.3 2 0 48 
76 0.45 : 0 92 


while such copper as is excreted by ananassae is lost quite quickly and the 
fraction retained is rather higher. The significance of this will be discussed 
further below. 

Although a variety of experiments were tried using food enriched with 
different amounts of stable copper for flushing out, no evidence was ob 
tained for an increased rate of excretion under these conditions. This is 
illustrated in lines 2 and 3 of table 1. The line marked “Cu” represents 
larvae scoured on food containing 25 wg. Cu per gram, whereas the other 
groups all were on plain food. 

To establish the distribution of copper among tissues and organs both 


counting and autoradiographic techniques were applied. For counting 
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larvae were dissected into fractions consisting of the middle midgut, the 
rest of the gut, and the remainder of the larva. In some instances the gut 
was further subdivided and the malpighian tubules were separated. Since 
considerable difficulty is entailed in recovering the malpighian tubules im 
foto they were usually included with the gut remnant instead of being 
counted separately. 

The histochemical investigations cited above had indicated that while 
there is a more or less constant quantity of detectable copper in malpighian 
tubule cells, most of the copper which can be revealed in copper-fed larvae 
1s located in the cells of a small section of the middle midgut. From the 
counting data set out in table 5, it may be seen that this predominance in 
histochemical copper was illusory. Doubtless, the sources of error in this 


FABLE 5 


Cu Uprake BY DROSOPHILA LARVAE AND [tS DISTRIBUTION IN RELATION ro Cu LEVEI 
CU IN uG. CU/LARVA % IN % IN % IN 

MEDIUM AFTER 8S HRS. ON MID-MID KEST OF REST OF 

SPECIES ; ; AND 4-8 HRS OFF Guy au? LARVA 


Repleta 2: 10 X 107 16 56 27 
105 & 10 4 3.é 15.7 41 
1406 «Kk 10 ? 38 1 21 
2260 xX 10 4 ‘ 37 

5 1074 1 § 
109 « 10 54 
$57 10 

H04 * 10! $5.§ 


Ananassac 


Melanogastetr 2! j Tie ‘ 17 
: 31 * 10 


762 10 & 59S 
$14 & 10 ) 


* Very little data 


case were both concentration (that only in this region was the copper con 
centrated enough to give a visible color reaction) and masking (that a 
larger part of the non-middle mid-gut copper is in combinations not reac 
tive in the test used). In the case of the rep/eta larvae the counting rates 
were enough above background and the statistical significance high enough 
to provide complete confidence in the reality of the increase in proportionate 
amount of middle mid-gut copper, as the copper supply in the medium was 
increased from 2.5 to 200 wg. per gram. In the cases of ananassae and 
melanogaster this increase is also apparent from 2.5 to 100 ug. Cu per gram. 
Further data are needed to establish the difference between ananassae and 
melanogaster at 200 ug. Cu per gram, in that even though both show 
toxic effects by reduction in Cu intake, in the former the middle mid-gut 
fraction has also been reduced whereas in the latter it has continued to in 


crease. 
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Although the counting rates were quite low, we are inclined to regard as 
real, the uniformly greater proportion of middle mid-gut copper im larvae 
on 0.25 ug. as opposed to 2.5 ug. per gram. 

In Figure | are illustrated representative whole-mount autoradiographs 
of dissected larval guts. In these the photographic exposure time was 
uniform so that differences in degree of blackening closely approximate 
differences in content of Cu®. It is evident that small amounts of copper 
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FIGURE 1 


Positive prints of contact autoradiographs showing Cu®! in gut epithelium and Mal- 
pighian tubules. Anterior ends above, */,) natural size 


are distributed throughout the dissections even in the distal portions of 
the caeca and in the malpighian tubules although the greatest concentra 
tion is clearly in the middle mid-gut. It is planned to present elsewhere 
an extended discussion of our combined histochemical and histological and 
contact autoradiographic studies. The data presented here, however, 
require to be supplemented with the statement that these other studies 
indicate: 
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(1) Copper excretion to be largely, if not exclusively, from the cells of 
the middle mid-gut. The excreted copper is not found to penetrate the 
peritrophic membrane but passes slowly back to the hind-gut within the 
space between the peritrophic membrane and gut epithelium. The be- 
havior of copper in this space below the middle mid-gut, as had that of 
barium above the middle mid-gut, indicates a slow posteriad flow of fluid 
between peritrophic membrane and mid-gut epithelium. This middle 
mid-gut excretion continues over a much longer span of time in larvae of 
repleta than in ananassae. 

(2) A larger proportion of the total copper to be located in the malpig 
hian tubules as the copper in the food approaches normal levels. 


PABLE 6 


Cu® UpraKE BY DROSOPHILA FROM RADIOACTIVE YEAST SUSPENSION. EFFECT OF 
STABLE Cu CONCENTRATION IN SUSPENDING MEDIUM 


cu CONC RANGE RANGE 
DROSOPHILA VEAST IN SALINE © MIN. /LARVA Cc MIN. /LARVA 
SPECIES GENUS ua./G AFTER 10 HRS AFTER 19 HRS 


Melanogaster ‘andida 0 117-578 90-186 
Melanogaster ‘andida 2.8 63-121 54-119 
Melanogaster ‘andida 25.0 285-3905 89-237 
Repleta “andida 0 365-840 141-482 
Repleta ‘andida 2. 226-692 148-953 
Repleta Candida 25 368-459 


Replet: Debaryomyces 20-68 137 
Replet: Debaryomyces f 37-139 7-99 

Replet: Debaryomyces 2 9-116 149 
Replet: Hansenula 2! 11-480 5-396 


Cu® In SALINE ONLY 


Replet: Debaryomyces 5.0 6-307 
Replet: Debaryomyces 25.0 151-1141 83-1398 


252 


(3) Some copper in the stored contents of the ascending enterior 
malpighian tubules, in repleta larvae. 

(4) The distribution of copper to tissues aside from the intestinal tract, 
to be extremely sparing at low copper uptake, and never very large. 

In another series of experiments we investigated the larval uptake of 


copper contained in yeast cells. Cu*! was added to the liquid media on 
which yeasts of the genera Candida, Iansenula and Debaryomyces were 
growing and after time for uptake, the yeast cells were centrifuged out, 


washed with distilled water and resuspended in saline solution. This 
suspension was placed on filter paper and the larvae allowed to feed on the 
cells. The data presented in table 6 indicate not only that’ there was 
rapid and considerable absorption of Cu®! from the bodies of ingested 
yeast, but that in some way this Cu® is absorbed without mixing with 
stable copper added to the saline in which the yeasts were suspended. 
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The lower lines of table 6 indicate that Cu®! in saline with non-radioactive 
yeasts is readily absorbed by the larvae. Table 2 has already presented 
data to show that at 25ug. Cu per gram, absorption is no longer propor 
tional to concentration, so that it does not appear possible that the lack of 
dilution effect shown in table 6 may be explained by uptake of all the 
copper passing through the digestive tract. 

Actually it appears that both chemical and morphological differentiation 
exists between yeast copper and medium (or perhaps soluble) copper. 
Histological autoradiographs of larvae at various times after feeding on 
Cu®! yeast show that the copper is not removed from the yeast cells until 
they have passed into the posterior mid-gut and that copper taken up from 
these is not concentrated in the middle mid-gut as ‘‘soluble’”’ copper is. 

By taking the proportionality of uptake to copper concentration in the 
food, established in table 2, correcting to the plateau for the time as in- 
dicated in table 3, and to zero excretion as indicated in table 4, we arrive at 
the following factors for copper content in late third instar larvae. These 
are expressed as wg. Cu per whole larva, or per mg. dry larva, per ug. Cu 
per gram of the medium on which the larvae have developed. 

The factors are 

For repleta 6 & 10°? wg. Cu per larva or 6 & 10 * per mg. dry wt 
For ananassae 7 X 10 4 yg. Cu per larva or 14 & 10 4 per mg. dry wt 
For melanogaster 5 & 10 3 yg. Cu per larva or 14 & 1074 per mg. dry wt 


The figure for melanogaster assumes that its uptake and excretion curves 
are quite like those of ananassae. Our analyses for the copper content of 
Drosophila media are quite variable but indicate a range in the neighbor 
hood of 0.50 ug. Cu per gram fresh weight. 

In the case of larvae raised on the two commonest fly food yeasts, 
Candida and Saccharomyces, it appears from other studies to be reported 
elsewhere, that most of the copper in the medium will be taken in as 
“soluble” copper, and will be metabolized as such. Thus we should expect 
20%) and perhaps much more of the total copper of normal larvae to be 
found in the cells of the “Cu band” of the middle mid-gut. It is gratifying 
to be able to report that these figures are well within the analytical range 
for copper when approached by neutron capture activation analysis, and 
that the arrangements are almost complete to perform the required analysis. 

In other experiments, not reported in detail here, we have found that 
Drosophila tunebris, a member of the subgenus Drosophila, as is D. repleta, 
has about the same coefficients for copper uptake. As discussed else 
where! we have postulated on histochemical grounds a profound difference 
in copper metabolism between the two subgenera Drosophila and Sopho 
phora (which includes D. ananassae and D. melanogaster). The radio 


copper studies reported here support this hypothesis, as indicated above, 
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by showing differences in uptake, rate of excretion, and distribution within 
the organism, between members of the two subgenera. We hope that by 
extension of such studies as these to more closely related species of the 
genus and to some of the inter-specific and inter-racial crosses, it may be 
possible to cast light on the gene mechanisms controlling copper metab- 
olism. 

Summary. The uptake of copper by larvae of four species of Drosophila, 
traced by Cu", is proportional to the copper concentration in the medium 
over the range 0.25-10 K ug. Cu/g. Above this level uptake falls off. 
Rates of excretion as well as distribution in the tissues have also been de- 
termined. From these data factors relating copper content of larvae to 
copper concentration in the medium have been calculated. Both counting 
and autographic methods demonstrate that a large fraction of the tissue 
copper of Drosophila is in a form not demonstrable with presently available 
histochemical techniques. Further support is lent to the hypothesis, 
previously advaneed,* of a profound difference in copper metabolism be 
tween two of the major subgenera (Drosophila and Sophophora) of this 
genus. 

It is shown that Cu! ingested as part of the yeast cell is absorbed without 
an opportunity of mixing with stable ionic copper simultaneously ingested. 
Thus there appear to be at least two pathways of copper uptake, one for 
ionic, the other for bound forms. 
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CONFORMAL RELATIVITS 
By RiIcHARD L. INGRAHAM 
INSTITUTE FOR ADVANCED STUDY, PRINCETON, NEW JERSEY 
Communicated by O. Zariski, August 2, 1952 


If one keeps invariant angle but drops invariant (four-dimensional 
length from an essentially Einsteinian description of the world, one gets 
Conformal Relativity. This theory, skeletonized below, will be published 
elsewhere. 

The basic idea is to let the invariance of the light cone alone (without 
the additional requirement that there exist an invariant length) define the 
class of preferred systems (i.e., those observers for whom the nature laws 
are “in their simplest form’). The field equations of the resulting theory 
thus exhibit imzvariance in form under the whole conformal group in local 
space-time, interpreted as coordinate transformation group. Accordingly, 
we define the Special Theory of Conformal Relativity to be the study of the 
ordinary conformal (“conformal flat’’) space C, with angle-defining form of 
signature (+++ —) | when the four-dimensional manifold is interpreted 
as space-time. It treats the kinematical relationship between equivalent 
observers and the behavior of measuring apparatus in the absence of all 
“force” fields. 

The mathematical techniques needed are classical, and some of the re 


sults are already in the literature. The main result is due to J. Haantjes' 


who in 1940 discovered that the conformal group represents transforma 
tions between all observers (cartesian systems in different Ry.) 1n uniform 
relative acceleration. All relative velocities are bounded above by light 
velocity, which is of course the same for all observers. For details on the 
non-invariance of length and rest mass,’ see Haantjes /oc. cit. Other results 
proper to the Special Theory are, e.g. the behavior of meter sticks and 
clocks under conformal transformations, the composition law of relative ac 

celerations, etc 
An event of this theory is defined to be a hypersphere (center x 
4), radius R), the subclass of nullspheres* (R = QO), identified with their 


m 


(m = l, 


center points x”, corresponding 1-1 with the events of ordinary relativity 


Hexaspherical coordinates,‘ the preferred systems of C,, refer these events 
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to six real independent hyperspheres as coordinate surfaces. Thereby 
points (nullspheres) of C, have hexaspherical coordinates Z” (a = 0, 1, 
5) satisfying a quadratic equation 


SnZ'Z = 0. (1) 


The subgroup of the group of all constant linear transtormations of the 


. P . cpa’ a ob ° 
hexaspherical coordinates of events Z° = P%,Z"(a’ = a) which takes 


points of Cy into points of Cy: 
Pr la P’, = hS,, (for some h) , (2) 


is just the 15 parameter group of conformal (coordinate) transformations 
of points of Cy." This is then the fundamental group of the Special Theory. 
The physical interpretation of the hypersphere-event is the specification 


~m 


of position x” in space-time and some standard reference length or gauge 
R (real or imaginary if space- or time-like, resp.) at that position. We 
have to do then, in Conformal Relativity, with a five-dimensional contin 
uum, the space-time-gauge continuum. (Gauge at a point in general changes 
under a conformal transformation, and the proportionate change is not the 
same for all gauges, unlike certain previously considered gauge theories.® 
Only zero gauge is preserved, 

Since the hexaspherical coordinates Z¢ of events in C, define an ordinary 
real five-dimensional projective space P;, everything can be equivalently 
phrased in the five-dimensional projective language in a well-known way.® 
I.e. (we give it for later convenience), the geometry of an ordinary confor 
mal space C,, (of signature Sig. g»,, Say) is identical with that of an ordinary 
real n + |-dimensional projective space P,,,; provided with a special quad 
ric S,» of signature (—1, Sig. gma, +1) representing points identified with 
nullspheres of X,, (ef. Ref. 1), whose group is restricted to those collinea 
tions taking the quadric into itself (ef. Ref. 2). 

The Special Theory cannot properly be called a relativity mechanics, 
for there are no laws of motion under various forces. These are given by 
the General Theory of Conformal Relativity. It 1s defined as the study of a 
general (“curvilinear’’) conformal space A, (to be defined infra) with angle- 
defining form of signature (+++ —) | over the X, of space-time-gauge 
which is subjected to the one global condition of minimum total (1.e., inte 
grated) curvature. In addition to the kinematical results known from the 
Special Theory, one gets as extremal equations a set of nature laws de 
scribing various ‘‘force’’ fields. 

We define: A general (‘curvilinear’) conformal space K, is ann + 1 
dimensional manifold X,,,, with ordinary conformal spaces C,, of the same 
signature as local spaces. In the n + 1-dimensional projective language 
(ef. the corresponding definition of C, supra), it is thus a curvilinear real 
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n + \-dimensional projective space //,,;" (homogeneous coordinates “ 
say) with special quadrie S,,(X“) homogeneous of some degree. 

For the purposes of the General Theory, we then start with a Ay whose 
manifold X, is interpreted as the space-time-gauge continuum and whose 
quadric S,, 1s of signature (—, +++-—, +)1. In addition it is by hy- 
pothesis linearly connected: we are given a projective connection [,,*. 
I," must be homogeneous of degree —1 so that covariant differentiation 
wrt it takes projective quantities of //,,, into projective quantities of 
/1,,,;.. Also in order that the total curvature (given by (3) infra) exist as a 
projective invariant, the degree of S,, must be chosen —2. Let the curva- 
ture tensor of I’ be N,,,-4, the Ricci tensor V,, == N,».", the Gaussian curva- 
ture .V S°N,». Then the condition of minimum total curvature re 
quires 


5S NV S(dX) = 0 (S = Det S,,, (dX) dX° dX"... .dX*) (3) 


As in Einsteinian relativity, this global condition is all that is required of 
the geometry. 

The connection I’ ard the quadric S are independently variable in (3). 
It turns out that I’ is necessarily not only symmetric but must be the Chris- 
toffel symbols of .S.° 
tation, Maxwell-type electromagnetism, and several mesons (tentatively inter- 
preted). Thus a unified field theory of force fields necessarily results from 


These extremal equations describe instein-type gravt- 


dropping invariant length from an essentially Einsteinian description of 
nature. ‘The mesons are vector and scalar wrt a group transforming four- 
position variables; they are particles of integral spin and non-zero rest 
mass (defined a posteriori of course*) in the language of particle theory. 
The electromagnetic and mesic theories thereby given are approximately 
linear in regions of tenuous energy, becoming essentially non-linear in 
regions where energy is dense. The interpretation (classical, quantum 

) is left open. The field equations are covariant against the group of 
all holonomie and non-holonomic transformations in //;, that 1s, under 
arbitrary projective coordinate transformations in the local P;’s. In the 
four-dimensional language, they express the vanishing of a hexaspherical 
tensor which is formed out of the curvature in the familiar way. The field 


equations are also all conformally invariant in form (-1.e., maintain their 


“simplest form’) under the whole conformal group or, kinematically, with 
respect to the class of all observers in uniform relative acceleration. 

Conformal Relativity’s relation to the former relativities, projective and 
afline (or ‘“‘general’’), with which it forms a generic class, will be treated in 
detail in the article to appear; here only a few of the more important re 
marks can be made. 

1. In many ways, projective relativity'’ (or the Kaluza theory'') re- 
veals itself as a provisional stage on the way to Conformal Relativity in 
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that features which were anomalous or arbitrary in the former are only ex- 
plained or uniquely fixed when we come to the latter; e.g.: 

(a) The point of view advanced in the present theory is that the appear- 
ance of projective geometry in unified field theory arises naturally only in 
handling the conformal geometry of four-dimensions in the framework 
of a five-dimensional projective tensor calculus. 

(6) The anomalous 5 * 5 quadric of the projective relativity ——Kaluza 
theory, when completed to the 6 6 quadric (1) of this theory, has the 
well-defined geometric role of defining the point-events among the class of 
hypersphere-events. Thus incidentally, any skew part of S,, would be 
geometrically meaningless here, and is therefore excluded from considera- 
tion,!” 

(c) There is no need for any arbitrary restriction on spatial symmetry of 
the cylinder condition type to eliminate an unwanted dimension-— space- 
tine-gauge fills out a whole Xs. 

(d) A priori criteria were lacking to define completely the quadric’s sig- 
nature in projective relativity.'* But it is a fact that the quadric’s signa- 
ture in this theory is fixed as (—1, Sig. gm,, +1) by the signature of space- 
time plus the demand that the local hexaspherical coordinates be real 
i.e., that the local spaces be real Ps's. 

2. Fundamental is the fact that only in this theory do we get the iden 
tity of the class of all observers “‘indistinguishable with respect to the light 
cone’ (kinematically, all those in uniform relative acceleration) and the 
class of all observers “‘indistinguishable with respect to the nature laws” 
(preferred systems). Calling these classes C; and C:, respectively, in Con- 
formal Relativity C; = C.. In projective relativity, C, is not a subset of 
Cs, Cy is not a subset of C). In affine (‘general’) relativity, C, ¢ Cj, but 
C, # C;. In both latter relativities, the intersection C; N C, is a class of 
observers connected by orthogonal transformations, so that an invariant 
length figures in both theories. 


' Haantjes, J., Proc. Ned. Akad. v. Wet., 43, 3 (1940) 

2 Haantjes treats only the classical Maxwell theory, where rest mass is of course rep 
resented by a number, a priori given, attached to each particle. 

‘In view of the signature (+++ —) 1, nullspheres are light cones, viewed as real 
loci, and centers x” are vertices 

‘For details on polyspherical coordinates, see Klein, F., Vorlesungen tiber Hohere 
Geonietrie, 3rd ed., Springer, Berlin, 1926 

®’ Weyl, H., Raum. Zeit. Materie, 4th ed., Springer, Berlin, 1921 

6 Klein, F., Math. Ann., §, 257 (1872) 

? For details on the curvilinear projective geometries /7,, see Schouten, J., and Haant 
jes, J., Comp. Math., 3, 1 (1936) 

* Cf. Schrodinger, E., Proc. Roy. Trish Acad., 51A, 63 (1947) for the four-dimensional 
affine case. 

® Cf, Einstein, Infeld, Hoffmann, Ann. Math., 39, 65(1938). And it is worth remark 


ing that the non-invariance of this rest mass under conformal transformations will be 
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Vor 


the explanation offered by Conformal Relativity of the apparent multiplicity of certain 


kinds of integral spin particles encountered experimentally 
" Veblen, O., and Hoffman, B., Phys. Rev., 36, 810 (1980) 
' Kaluza, T., Sitsher. Preuss. Akad. Wiss., 966 (1921 
'2 In contrast to Einstein, A.,*/bid., 414 (1925), or (8) 
'§ Pauli, W., Ann. Physik, V18, 305 (1933) 
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